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We consider a Markov evolution of lozenge tilings of a quarter-plane 
I— I and study its asymptotics at large times. One of the boundary rays serves 

^-H as a reflecting wall. 

We observe frozen and liquid regions, prove convergence of the local 
• correlations to translation-invariant Gibbs measures in the liquid region, 

and obtain new discrete Jacobi and symmetric Pearcey determinantal 
d point processes near the wall. 

Ch The model can be viewed as the one-parameter family of Plancherel 

I I measures for the infinite-dimensional orthogonal group, and we use this 

interpretation to derive the determinantal formula for the correlation func- 
tions at any finite time moment. 
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1 Introduction 

The principal object of study in this paper is a one-parameter family of proba- 
bility measures on certain interlacing two-dimensional particle systems that can 
be defined in at least three different ways. 

Random lozenge tilings. Consider the domain pictured on the left in Figure 
[ijdrawn on the regular triangular lattice, and consider all possible tilings of this 
domain by lozenge^ 



Figure 1: Lozenge tiling and nonintersecting paths 




A lozenge consists of two neighboring elementary triangles glued together. 
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An example of lozenge tiling can be seen in the middle of Figure [T| To each 
tiling we assign a weight equal to ^ raised to the number of vertical lozenges 
on the left border of the domain (three such lozenges are highlighted on the 
figure). Let us normalize the weights so that the total weight of all tilings is 
1; then we obtain a probability distribution with three parameters a, 6, c that 
represent side lengths of our domain. 

Let us further consider the limit a, 6, c ^ oo so that 2ab/c — > t > 0, and focus 
on the part of the tiling that is of finite distance to the bottom-left corner of the 
domain. One can show that in this limit our probability distributions weakly 
converge to a probability measure M.t on lozenge tilings of the quarter-plane, 
and it is the limiting measure that we are interested in. 

Lozenge tilings are also commonly viewed as stepped surfaces (when three 
types of lozenges are interpreted as three faces of 1 x 1 x 1 cubes in a three- 
dimensional space) , as nonintersecting paths (see the right-most part of Figure 
[T]) , and as dimers on the hexagonal lattice (see Figure |4] in Section |2.3| below) . 
Theory of dimer models is a rapidly developing subject, see |29j for a recent 
review and references. 

In terms of nonintersecting paths, the initial (a, b, c)-measures give an extra 
factor of 2 every time the left-most path passes by the wall. Thus, it is natural 
to say that this path reflects off the wall. 

Random surface growth. Any lozenge tiling is uniquely determined by loca- 
tions of lozenges of a single type. Let us introduce coordinates on the plane as 
shown in Figures [T] and |4] and mark the midpoints of all vertical lozenges; call 
them particles. Denote the horizontal coordinates of all particles with vertical 
coordinate m by y™ > 2/™ > . . . . Then f is a probability measure on particle 
configurations 

|fc=l,2,...,[^]; m = l,2,...}cZ>o 

that satisfy the interlacing conditions y^^^i < y™ < y™^"'^ for all meaningful 
values of k and m. 

We show that Ait is the time t distribution of a continuous time Markov 
chain defined as follows. 

The initial condition is a single particle configuration when all the particles 
are as much to the left as possible, i.e. y™ = m — 2A: -I- 1 for all k, m. Now let 
us describe the evolution. 

We say that a particle y™ is blocked on the right if y™ + 1 = y'^Zi , and it 
is blocked on the left if y™ — 1 = y]^~^ (if the corresponding particle y^Si or 
y™~^ does not exist, then y^ is not blocked). 

Each particle has two exponential clocks of rate ^ ; all clocks are independent. 
One clock is responsible for the right jumps, while the other is responsible for 
the left jumps. When the clock rings, the particle tries to jump by 1 in the 
corresponding direction. If the particle is blocked, then it stays still. If the 
particle is against the wall (i.e. yTL+ii = 0) and the left jump clock rings, the 
particle is reflected, and it tries to jump to the right instead. 
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When J/™ tries to jump to the right (and it is not blocked on the right), we 
find the largest r G Z>o U {+00} such that y™^* = + i for < i < r, and 
the jump consists of all particles {?/r^'}i-o moving to the right by 1. Similarly, 
when y™ tries to jump to the left (not being blocked on the left), we find the 
largest I G Z>o U {+00} such that J/j^"' = y]^ — j for < j < Z, and the jump 

consists of all particles {y^'^i^'j"' moving to the left by 1. 

In other words, the particles with smaller upper indices can be thought of 
as heavier than those with larger upper indices, and the heavier particles block 
and push the lighter ones so that the interlacing conditions are preserved. 



Figure 2: First three jumps 




Figure [2] depicts three possible first jumps: Left clock of yl rings first (it gets 
reflected by the wall), then right clock of yf rings, and then left clock of yl 
again. 

In terms of the underlying stepped surface, the evolution can be described 
by saying that we add possible "sticks" with base 1x1 and arbitrary length of 
a fixed orientation with rate 1/2, remove possible "sticks" with base 1x1 and 
a different orientation with rate 1/2, and the rate of removing sticks that touch 
the left border is doubled 



Figure 3: Adding and removing "sticks' 



Add 




^This phrase is based on the convention that \iv is a figure of a 1 X 1 X 1 cube. If one 
uses the dual convention that this is a cube-shaped hole then the orientations of the sticks to 
be added and removed have to be interchanged, and the tiling representations of the sticks 
change as well. 
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A computer simulation of this dynamics can be found at 



http : //www.math. caltech.edu/papers/Orth_Planch. html 



Similar Markov chains have been previously studied in [8J without the wall, 
and in [39] with a different ( "symplectic" ) interaction with the wall. 

Representation Theory. Let 0{N) be the group of N x N orthogonal matri- 
ces with real entries. The group 0{N) is embedded in 0{N + 1) as a subgroup 
of matrices fixing the {N + l)st basis vector. Let 0(oo) = lJw'=i ^i-^) 
infinite-dimensional orthogonal group. 

The measures Mt are the Fourier transforms of the distinguished one-para- 
meter family of indecomposable characters of 0(oo) (the indecomposable char- 
acters of 0(oo) were classified in '38] as a part of a solution of a much more 
general problem). It is natural to call them the Plancherel measures. Details 
can be found in Section [21 

Similarly defined Plancherel measures for the infinite symmetric group S{oo) 
and the infinite-dimensional unitary group U{oo) have been thoroughly studied, 
see [SgilllllHllalllllMlinilMlEa] for 5(oo) and [SlllSllII] for U{oo). 

Results. We first prove, see Theorem |3 . 1 2| below , that representation theoretic 
and Markov chain descriptions of M.t given above are equivalent (the lozenge 
tiling description oi M.t is a simple corollary of the representation theoretic one 
and Theorem 1.4 of |38jV This equivalence is far from being obvious, and we 
employ the general formalism of ISl to give a proof. 



Our second result (Theorem 4.11 shows that M.t, viewed as a measure on 
particle configurations {y^}, is a determinantal random point process (see Ap- 
pendix |A] for basic definitions), and it also provides an explicit formula for the 
correlation kernel. In fact, we prove such a result for random point processes 
associated with arbitrary indecomposable characters of 0(oo). 

We then focus on the asymptotics of A^t as i ^ oo. 

As one might anticipate from previous results on dimer models and Plancherel 
measures, cf. [281 ISHl HI E], as t — > oo the quarter-plane should spHt into 
"frozen" parts and a "liquid" part. In each frozen part the tiling asymptoti- 
cally consists of lozenges of only one type, while in the liquid part the random 
tiling locally (i.e. on the lattice scale) converges to the unique (thanks to [3^ 1 
translation invariant Gibbs measure of a certain slope; the slope depends on 
the location in the liquid region. The underlying random surface should also 
converge, in a suitable metric, to the deterministic smooth limit surface, and 
the slopes of the Gibbs measures are the slopes of the tangent planes to this 
limit shape. 

In Theorem 5.2 we prove the statements about local convergence. The frozen 
and liquid phases can be clearly seen on Figure [S] More exactly, we prove the 
convergence of our correlation kernel to the incomplete beta-kernel first obtained 
in [27J [39], see [T3] for a detailed discussion of the Gibbs properties of the 
corresponding determinantal process. In Section [5. 2 [ we also provide a formula 
for the hypothetical limit shape, although we do not address the concentration 
of measure phenomenon. 

From previously known results it is also natural to expect that near the 
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boundaries between frozen and liquid regions away from the wall, our determi- 
nantal process converges in an appropriate scaling to the so-called Airy process, 
see e.g. Section 4.5 of jTT] for an analogous results in the case of Plancherel 
measures for U{oo). This is indeed correct, and since the result and the method 
of proving it are well known by now, we did not include them in this paper. 

The main novel feature of the model analyzed in this paper is the wall, and 
we focus on the corresponding scaling limits. 

The simplest case is the neighborhood of the origin (the corner of the quarter- 
plane) . Taking t ^ oo asymptotics in Theorem |2.8[ one can easily show (al- 
though we do not do this in the paper) that as one scales the horizontal co- 
ordinate by and keeps the vertical coordinate finite, A4t converges to the 
antisymmetric GUE minor process (aGUEM) of see also [TSl [1^]. Note 
that the way this process was obtained in [22] from lozenge tilings of a half- 
hexagon is also similar to what we are doing. The aGUEM process can also be 
obtained from the evolution of interacting Brownian motions with a reflecting 
wall, which can be seen as a limit of the Markov chain described above; see 
[Smn] for details. 

The first genuinely new limit that we obtain takes place in the region where 
the liquid part meets the wall. In Theorem |5.6| we show that on the lattice 
scale our determinantal process converges to a limiting determinantal process 
on Z>o X Z that is translation invariant in the second coordinate. We use the 
term discrete Jacobi kernel for the correlation kernel of the limiting process. 

The second new determinantal process arises when we look near the loca- 
tion whe re th e boundary between frozen and liquid phases meets the wall. In 
Theorem 
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we prove that as one scales the horizontal coordinate by and 
the vertical one hy t^ , the correlation functions of our point process converge 
to the determinants of the kernel /C((Ti, cr2, 772) on M+ x M defines as follows. 

Let C be the contour in C consisting of rays from ooe^^^^ to to ooe~^'"^^. 
Then 



/C(ai,77i;a2,772) = ^ / / e-''^^'+''^"'-^'+"%os(aix) cos(a27.) 



TT^Z J J — 

2v/vr(r7i - r,,) Mm - m) ^ 4(,y2 - .71) ) 

We call it the symmetric Pearcey kernel because of the similarity of the 
above expression to the Pearcey kernel that has previously appeared in [21 
[m [m HOJ [46] . Using the nonintersecting paths interpretation mentioned above, 
it seems plausible that the symmetric Pearcey kernel should also appear in the 
model treated in [53] at the critical location when the paths touch the wall, but 
so far we have not pursued this connection any further. 

Acknowledgements. The authors are very grateful to Grigori Olshanski 
for a number of valuable remarks. The first named author (A. B.) was partially 
supported by the NSF grant DMS-0707163. 
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2 Measures on partitions 



2.1 Representations of Orthogonal Groups 

Let 0{N) denote the group of all real-valued N x N orthogonal matrices. For 
each N, 0{N) is naturally embedded in 0{N + 1) as the subgroup fixing the 
{N + l)-st basis vector. Equivalently, O G 0{N) can be thought of as an 
(iV + 1) X (A^ + 1) matrix by setting Oi,Ar+i = On+ij — for I < i,j < N and 
On+i,n+i = 1- The union lJw'=i ^i^) denoted by 0(oo). 

Let us review some basic results from the representation theory of finite- and 
infinite-dimensional orthogonal groups, see e.g. [38] . 

A character of 0(oo) is a positive definite function x '■ 0{oo) C which is 
constant on conjugacy classes and normalized, i.e. x(e) = 1. We further assume 
that X is continuous on each 0{N) C 0(oo). The set of all characters of 0{oo) 
is convex, and the extreme points of this set are called extreme characters. 

The set of extreme characters can be parametrized. Let M.°° denote the 
product of countably many copies of K. Let O be the set of all (a, /3, S) such 
that 

a = (ai > as > . . . > 0) e ]R°°, /3 = > /^z > ■ • • > 0) e 5 gR, 

OO 

5^(az + A) < 6. 

i=l 

Set 

OO 

1=1 

The special orthogonal group, denoted by SO{N), is the subgroup of 0{N) 
consisting of matrices with determinant 1. Let O E SO{N). li N = 2m is 
even, then the spectrum of any O is of the form {zi, z^^ , . . . , Zm, z^^}, while 
if A'' = 2m -I- 1 is odd, then the spectrum of any O G SO{N) is of the form 
{zi, zj~^, . . . , Zm, z~^, 1}, where in both cases Zi are complex numbers having 
absolute value 1. For this paper, if x is a character of 0{oo), SO{2m) or 
SO{2m + 1), then x{0) is written interchangably with xi^ii ■ ■ ■ t ^m)- Using 
this notation, any w G defines a function on 50(oo) = IJn'=i S0{N) by (see 
Theorem 1.4 of f31]) 

where 

^ ,5 (.+.--2) n (i + |(z-i))(i-t-§(z-i-i)) 

2 ; 
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or by setting x = {z + z ^) /2, 



■yi.-i) l-A(l-x) + /3f(l~a:)/2 
+ a^(l-x) + a2(l_a;)/2' 



i=l 



Note that the infinite product converges because J^i'^i + Pi) finite. As uj 
ranges over f2, the functions x'^ range over all the extreme characters of 0(oo) 
(Theorem 5.2 of [35]) • 

A partition of length < N is a, sequence of nonincreasing nonnegative in- 
tegers A — (Ai > . . . > Xn > 0). It is a classical result that the set of all 
irreducible representations of SO{2N + 1) over C is parameterized by partitions 
of length < N. The character of the irreducible representation of S0{2N + 
1) parameterized by A will be denoted by Xso{2N+i)' dimension by 

dim5o(2Af+i) X- Similarly, the set of all irreducible representations of SO{2N) 
over C is parameterized by sequences of integers A = (Ai,... Xn) satisfying 
Ai > . . . > Aat-i > IAatI, and the corresponding character and dimension are 
denoted by Xso{2N) ^^"^ dim5C)(2Af) A. Let In denote the set of all partitions 
of length < (J stands for Jacobi, see below). For convenience of notation, let 
A* = (Ai, A2, . . . , Xn~i, -Ajv) for any A e J^v- 

For any a; G fi, the restriction of x'^ to any SO{M) defines two measures 
Pn.-i/2 Pn^i/2 on by 

^uj, _ puj Xso{2N+l) , . 

X \sOi2N+l) - 2^ Pna/2[^) 1 ' W 

Ate dimso(2Ar+i) A 

^c., _ Sr- pu^ l^^ Xso(2N) + Xsoi2N) , . 

- .tC 2dim,o(2A.)A ' 

where x'^ is the character of 50(2A^+ 1) or SO{2N) parameterized by A. Evalu- 
ating both sides of the equation at the identity of the group shows that the sum 
of the weights is one. Furthermore, the weights are nonnegative because the 
characters are positive definite, so we obtain probability measures. Note that 
if the parameters of u are a = 0, f3 — 0,j ^ 0, then is the delta measure 
supported at the partition (0, 0, . . . , 0). 

There is a useful explicit formula for x^- Let j\^'^\x) denote the fc-th Jacobi 
polynomial with parameters a, 6 see e.g. |43J. Define the constant cu to be 

if ^>o, 

Cfc = < 2 • 4 • . . . • 2fc ' 

1, if k = {). 



and let i'-^'^Xx) = j\^-^\x)/ck. The character of O e SO{2N) or SO{2N+l) 
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IS 



det 



Xso(2Ar+l)(^ii ■ • ■ ' ^n) 



,(1/2-1/2) ( Zj+z-^ 



N 



det 
det 



N 



(XsO(2JV) + XsO(2Ar))(^li • ■ • i ^n) 



2J 



(-1/2-1/2) f 
\i-i+N 



N 



det 



_1N JV-i' 



Expressions Q and ([s]) can be simplified using 

,(-1/2-1/2) 
-■s 

j(l/2 -1/2) 



Z + Z 

z + z^ 



z + z 
2 ' 

^s+l/2 _ ^- 



s-1/2 



1/2 _ r-1/2 



The Jacobi polynomials also satisfy 



(a -1/2) 
k 



(^) = 2^1+7 



1 ,(a,-l/2) , 1 |(a,-l/2) 



2^-1 



±-, A;>0 



(1/2.-1/2) 



(1/2.-1/2) 



i^) + 2^1 



(1/2.-1/2) 



(-1/2.-1/2) 



(X) = Ji 



(-1/2,-1/2) 



Explicit formulas for the dimensions are 



dim5o(2Ar+i) A — 



n 

l<i<j<Af 



;2 - /2 

* J l<i<Af 



n 



where h — Xi + N~i+^, rrii — N ^ i 



and 



dim5o(2Af) A 



n 

l<i<j"<Af 



9 9 5 

^ J 



where li — Xi + N ~ i, nii — N 



Let /i^"'''^ denote the squared norm of 



MM) 



4'''''\xf{i-xr{i+x)''dx. 



Then 



hi^'"^ = 7rcl/W^^''\k), for a = ±i, 6 = -i. 
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where 



2, if /c > 0,a = = -i 



W('''''\k)= {l, if k = 0,a = b=-^ 

1, if k> 0,a= i,6= -i 
For proofs of these equations, see §1 of [37J, Chapter 4 of [43 , and Chapter 24 

of m\. 



In Section 2.5 a formula for Pj^ ^ will be proved. 



2.2 Central Measures 

For a — ±1/2, the measure on Jjv can be extended to a more general 

measure P". The purpose of this section is to explain how is constructed. 

Let In - and In,+ be two copies of In- Set J = UAr>i('If-'v,- U In,+)- Turn 
J into a graph as follows. Draw an edge between A G In,- and fi e Jjv,+ 
if < Xn < Mjv < • • ■ < -^1 < Ml- Draw an edge between A g In,+ and 
/X e In+1,- if /iw+i < -^w < ^J'N < Aat-i < . . . < Ai < /ii. It will be convenient 
to set Aat+i = 0, which gives the additional inequality Ajv+i < fJ-N+i- In either 
case, use the notation X ~< fi. 

Note that A ^ is equivalent to the following relation from representation 
theory. Let be the representation of SO{M) corresponding to /i and let V\ 
be the representation of SO{M — 1) corresponding to A. With this notation, V\ 
is a subrepresentation of V^lso(M-i) iff A ^ /x. See [50] . 

For any (A, /i) G JAr-i.+ x J^v,-, set 



>f(A,/x) 

If (A, ^) e Jat^- X In,+ , set 



1, if A ^, 
0, if A 7^ n. 



2, if A ^ ^, Ajv > 0, 
>c{X,^i) ^ {l, if X ^ Ajv = 0, 
0, if A 7« ^. 

The definition of x{X, fi) is motivated by the branching rules 

{Vx © no © 14, A e In,-, fi e Jw.+ (12) 

Ajv>0 Aiv=0 

^^150(2^-1) = ^Vx, X e In-i,+, fJ- e In,-- (13) 

A^/j 

A path in J is a sequence t = {t^'~ ^ t^'~^ ~< t^'~ -<...) such that S Jj,_ 
and € Ji_+ for 1 < iV. Let Jpat/ts denote the set of all paths in J. There are 
also finite paths, which are sequences u = {u^'^ -< -< u^'^ -<...-< u^'^) 
that end at some u^^*. Given such a finite path u, define the cylindrical set 
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Also for a finite path, define the weight Wu to be 

i/7„ = U^'^)>c{u^''^ , .... 

The brancing rules imply that if we sum Wu over all finite paths u that end at 
A, we get dim A = dim Vx- Note that dim V\ = dim Vx* . 
A probability measure | • | on Ipaths is called central if 

\Cu\ ^ lai 

Wu Wy 

for any two finite paths u, v that end at the same partition. For a more general 
definition of central measures, see section 6 of [32 . 

Each Lj E defines a central measure P'^ on Ipaths as follows. The measure 
of any cylindrical set Cu is given by Pjsj ±{um,±) / dimuAr-i-. This measure is 
well-defined because of the consistency relations (cf. ([2])-([3])) 

^ Pj^^+(/i)x(A,Ai)-Pj^,_(A), 

/iGJiv,+ 

2.3 Equivalent Interpretations of Paths in J 

It will be useful to interpret measures on Ipaths as random point processes on a 
two-dimensional lattice or as random lozenge tilings of a quarter plane. 
Set 

X = Z>o X Z>o X {±i}, 2} = Z>oxZ>o. 

For {nj,aj) G Z>o x {±|}, j — 0, 1, we write {uq, gq) < (ni, ai) if 2no + ao < 
2ni + ai. Also write (no, oq) < (ni, ai) if 2no + uq < 2ni + ai. That is, 

(l,-l/2)<(l,l/2)<j(2,-l/2)<.... 

Furthermore, set d(no, oq; n-i, ai) = |2(rii — no) + ai — ao\. In other words, 
(i(no, ao; ni, ai) is the distance between the levels (no,ao) and (ni,ai). 
We identify X and 2) via the bijection 

t:X^2), {x,n,a) 1-^ (2x + a + ^,2n + a — ^) . 

To any finite or infinite path A = (A'^-*'"^/^ -< A*-^''^/^ in J, we associate 

two point configurations (subsets) £-xW C X and >C2)(A) C 2) as follows: 

£3c(A) = {(4"^'",n,a) : 1 <A;<n, a S {±i}, n > l} , 
£2)(A) = {(yf+""'/',2n + a-i) :l<fc<n, a e {±i}, n > l} , 
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where 

^(„),a ^ ^(„),a ^ ^ _ ^.^ ^2„+a-l/2 ^ 2(a(")''^ + ,^ - fc) + a + i. 

Note that l{Cx{^)) — 'C2)(A) for any A. 

In Figure |4] black dots mark the elements of for A — ((1) -< (1) -< 

(1,0) (2,1) ^(2,1,1) ^(3,1,1)). 



Figure 4: 




The interlacing property for paths in J turns into 

This construction gives rise to a bijection between infinite paths in J and 
certain lozenge tilings of the quarter plane with boundary as indicated in Figure 
|4] Elements of CrQ{\) correspond to centers of lozenges of one specific type, as 
in Figure |4] The image oi^ipaths consists of lozenge tilings with [^^^^] lozenges 
of this type on the mth horizontal row, for m > 1. 

Equivalently, lozenge tilings can be viewed as dimers on the dual hexagonal 
lattice, see e.g. [IH]- In this language, paths in J correspond to dimers with 
exactly [^^^] vertical edges crossing the mth horizontal line. 

Thus, measures on ^ paths yield random point processes on X and 2), a mea- 
sure on lozenge tilings, and a measure on dimers. The centrality of such a 
measure on dimers can now be phrased in the following Gibbs-like manner: As- 
sign a weight of 1/2 to all the vertical edges of the hexagonal lattice that lie on 
the veritcal line with coordinate 0. These are marked by dashed lines in Figure 
|4[ All other edges have a weight of 1. Then, given the set of [^^^^] vertical 
edges crossing the mth horizontal line, the conditional distribution of the dimers 
below this line is proportional to the product of the edge weights of the dimers. 



Recall from sections 2.1 and 2.2 that any w G defines a probability measure 
P"^ on Ipaths- Let and denote the resulting point processes on X and 
2), respectively. If the parameters of uj are such that all ai — Pi = 0, then let 
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p1 y_ and p1 denote the fcth correlation functions of T"^ and , respectively. 
See Appendix |A] for general definitions of point processes. We will need the 
correlation functions in section [H 



2.4 Preliminary Lemmas 

Before continuing, a couple of lemmas will be needed. Since they will be used 
several times throughout this paper, it is convenient to gather them in this 
section. Lemma [2. 1| is a variant of the Cauchy-Binct formula. 



Lemma 2.1. For each nonnegative integer k, let and gk he some functions 
on C. For \ < i < N , let Xi, Ci be complex numbers such that X^fcLo ^k{Cj)9k{xi) 
converges absolutely for all I < i, j < N . Then 



J2 det[efc, (0)15.1 det[5fe, (x,)]5=i = det 

fci>...>/cjv>0 



.fe=0 



N 



Proof. The proof is almost identical to the proof of Theorem 1.2.1 of [25]. □ 

The next lemma is also useful. 
Lemma 2.2. For a = ±1/2, -1 < C < 1, and a test function T e C^[~l, I], 



E 

fc=0 



' J, 



w-^fc _ ^)a(^ _^ x)-'/^dx = r(c). 



(a,-l/2) 



Proof. Plugging in a = —1/2 and setting x = cos 0, C, = cos 9, (4.1.7) of [43] 
and (111 give 

^(-1/2,-1/2) 



jr 



(-1/2,-1/2) 



i if fc = 0, 

- COS k4> cos kO if fc > 0. 



Since T is C^, the Fourier series of T converges to T (see Chapter 3, Section 6 
of [S]): 

r(cos 0) To + f 1 cos + fa cos 2<?!) + . . . , 



where 



T{cos4>)d(l>, fc = 0, 

T{cos (j)) cos k(l)d(j), fc > 0. 
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Therefore 



OO ^1 



1 

TT 



2 r'^ 

T(cos0)c?0H — > / T (cofi (t>) cos k<j) cos kOdcf) 



To + Ti cos e + T2 cos 26 - 

no- 



In the case when a =1/2, (4.1.8) of f43^ and (|TT| tell us 

r(i/2,-i/2)|-^^^a/2,-i/2)(^) sin(fc + l/2)(j3 sin(fc + 1/2)61 



J, 



(1/2.-1/2) 



sm0/2 sin 6*72 



and the rest of the argument is similar. 

The previous two lemmas also imply the next one. 



□ 



Lemma 2.3. Let —1 < (i,...,(n < 1, and suppose ip G C^[— 1,1]. Let 
li, ... ,1m be nonnegative integers. Set 



where a = ±1/2. Then 



-1/2) 



AT 



J, 



(a -1/2) 



(x) 



j^K-l/2)ji/2 



^ det e^-'^'^iCi) det[5(fc„?.)]f,=i = ¥'(Ci)---^(C^)det 



-1/2) 



(C.) 



fei>...>fejv>0 

where 



N 



ij = l 



5(fc,0= / ekix)eiix)ip{x){l-x)''il + x)-^/^dx. 



Proof. By Lemma 2.1 



^ det [eg-^/^'CC.)]'^ det[5(fc„/.)]5=i = det 



fci>...>feiv>0 

By Lemma 12^21 



.k=0 



{Q9{k,h) 



E4"-^/'HC.)<?(fc,^,) = E / ei--^/^)(.)et-/^)(C.)<-^/^Ha:)^(.)(l-x)"(l + .)-V2d, 

it— n fc— n-'- 1 

(O)V'(C:). 



k=0 



k=0 ' 



^ g(i-l/2)^ 



□ 
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Lemma 2.4. The normalized Jacobi polynomials ji^^^^' ^^"^^ satisfy the follow- 
ing properties: 

(a) ^ M/(-l/2-l/2)(^)j(-l/2,-l/2)(^) ^ j(l/2,-l/2)(^)^ 



r=0 

s-1 .(-1/2,-1/2) 



(b) ^j(iA-i/2)(^) = ^ 



(X) 1 



X- 1 

(c) 1 /' j(i/2,-i/2)(2,)(i _ + = 1. 

J-i 

Proof. (a),(b) Let z be on the unit circle such that (z + z^^)/2 = x. Using (|6| 
and ([7|, the sum becomes a geometric series, which can be evaluated explicitly, 
(c) By part (a), the integral equals 

^ ^(-l/2,-l/2)(,) rl j(-l/2,-l/2)(^)j(-l/2.1/2)(^)(^_^)-l/2(^^^)-l/2^^^ 

which equals 1 by the orthogonality relations. □ 

Lemma 2.5. Let T G C^[— 1, 1]. The following identities hold: 
(a) 

g V^(-l/2,-l/2)(^) .1 _,(_i/2,-i/2)(^)y(^)(^ _ ^)-i/2(^ ^ = r(l). 

g ^(-l/2.-l/2)(^) r j(-V2,-l/2)(^)J.(^)(l„^)-l/2(i + ^)-l/2^^ 

- C }^'^^-^'^\x){T{l) -T{x)){l~ x)-^'\l + x)-^'^dx. 



r— s+l 



Proo/. (a) Note that 



1/2, 1/2) .1 j(-l/2,-l/2)(^)y(^)(^_^)-l/2(^^^)-l/2^^ 

oo 1 (-1/2,-1/2), . ^(-1/2,-1/2).-.. 

: E / . ^^^^ T^Tj^, — - + .r^'^dx 



-1 /ir 
=T(1) 

The first equality follows from (111 and JJ. (1) = 1, and the second 



equality follows from Lemma |2. 2 
(b) By Lemma [Q^c) , 



r(l) = - /' r(l)j(i/2-i/2)(a;)(i _ a;)-i/2(i + a;)-i/2d^. 
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By Lemma 2.4 'a), 



E 

r=0 



Subtracting this from the sum in (a) proves (b). 



□ 



Since the series in Lemma 2.5 a) converges, 
Corollary 2.6. For T eC^[-l,l], 



Hm 



j(-i/2'-i/2)(a;)T(a;)(l - x)-'/^{l + x)-^'^dx = 0. 



Lemma 2.7. For T e C^[-l, I], 

g V^-(l/2,-l/2)(^) rl j(l/2,-l/2)(^)y(^)(^_^)l/2(^^^)-l/2^^ 

1 o 1 



Proof. By Lemma 2.4 'b), 



^ ^(1/2, l/2)(^) .1 j(i/2,-l/2)(^)y(^)(i _ ^)l/2(i ^ :,)-l/2d^ 
„ — n ^ J — 1 



1 /"^ j(-l/2-l/2)^ 



Taking s — s- cx) and using Corollary |2.6| 
g Ty(i/^'-i/^)(r) j(i/2,-i/2)(^)y(^)(i _ ^)i/2(i ^ 

Q — n ^ -/—I 

1 /•! 

r(:E)(l-a;)-i/'(l + x)i/2da;. 



Subtracting these two sums proves the lemma. 

2.5 Explicit Formula for the Measures 

In this section, we prove the following statement: 



□ 



16 



Theorem 2.8. For any X = (Ai > . . . > Xn) G Jw, a = ±1/2; and lj efl, 

PN.aW = ^N.a ' dct [/j^'"^ (Aj - i + N)]i<,^j<N ' diniso(2JV+l/2+a) A. 



where 

AN. 



and 



(14) 



CN.a 



'2{N-i)N/2^ 1/ a= 1/2, 

2(^-2)(^-i)/2, z/a=-l/2. 



Theorem 2.8 follows from ([TJ and the following statement with E = . It 
is an orthogonal group analog of Lemma 6.5 of [4T]. 

Lemma 2.9. Let E{x) G C^[— 1,1]. Fix complex numbers zi,...,zjv on the 

unit circle and set (k ~ (zk + z'f^^)/2. Then 



i?(Ci)...£;(Civ) = 2.4.....2^-i 

X ^det[/f'^/^)(A.-z + A^) 
AeJiv 

and 

EiCi)-.-E{CN)^2-4- . 



l<iJ<N 



XsO{2N+l)i^l^ • ■ • ' ^n) 



J2det [/f^-^/^'(A.-z + iV) 



l<i,j<N 



■ (Xs'0(2Ar) + XS0(2N) ■ ■ iZn), 



where fj'^'°''{k) is defined by (14 1, with E{x) in place of E^{x). 



Proof. Let Ei{x) ~ x ^E{x) ior 1 <i < N. Using equation (111 and Lemma 
^2] shows that for 1 < i < iV, 

I — n-^ — 1 h. 



k=0 

oo 

E 

oo 



(NM)fl^^ l(a,-l/2) 
k 



k=0 

By Lemma [2.1 1 



<^-^/B(t)ji,"'"'/''(i)ji"'~'/')(a;)(l - tni + t)-'/^dt 
(x). 



^ det[/f''')(A,:-* + 7V) 



Ai>...>Aiv>0 



l<iJ<N 



■ det 



j(a -1/2) 



(0) 



l<i,j<N 



, (15) 
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where Xi ^ ki + i — N . Equation ( 15 1 can be rewritten as 



det [/f■'^)(A,-^ + iV) 



det 



.(a -l/2).> ^ 
-'Xi-i+N 



l<iJ<N 



(16) 



Set — (-Zfe+^^fe ^)/2. First consider the case when a — 1/2. Using Q, equation 
(16) becomes 



iJ(Ci)...£;(Civ) = 2-4.....2^-i 

X Edet[/f'^/^^(A.-^ + iV) 



Xso(2W+l)(^l; ■ • ■ j^Af)- 



Now consider the case when a = —1/2. Using (|5|, equation (16) becomes 



^det [/f-^/^'(A.-z + iV) 



2N) +XS0(2Ar))(^l) ■ • ■ ; ^w)- 



AeJiv 



l<i,j<N 



□ 



3 Stochastic Dynamics 

3.1 Markov Chain on One Level 

Let E{x) £ Ci[-l,l] such that £'(1) ^ 0. Plugging zi = . . . = Zn = 1 into 



Lemma 2.9 shows that E{x) defines a (possibly signed) normalized (i.e. all the 

dim.s 



weights add up to one) measure Pat.q on In. a by 



^V,a(A) = const • det 4^'°' (A, - i + N) 



l<i,j<N 



^SO{2N+l/2+a) 



A 



where /j^'"'' is defined by (|14|), with E instead of E'^ 



Fix <p G Ci[-l,l]. Define 



£:(a;) = E{x)ip{x) 



and define /j^'"' and P/v_a as in Theorem 2.8 except with E instead of E'^ 



For a = ±1/2, let be defined on Z>o x Z>o by 

oo 



k=0 



(a. ^^^-'(a;) and integrate over [—1,1]. Then the above 



Multiply both sides by J; 
definition is equivalent to 
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Proposition 3.1. With the above notation, 

PjV,a(A) 

dimso(2Af+i/2+a) A ^ — 



J2 det[I^i^i,-i+N, A,-j+7V)]i<,j< 



dini50(2Af+l/2+a) 

Proof. Let A, B, C be the following matrices: 

^<j<N, l<j<oo 
A{i,j)=I^{j,i), l<i,j<oo 
^(*, j) 1 < J < A^, 1 < « < oo 

Then C = AB, so by the Cauchy-Binet formula, 

det[/j'^'°^(A,-?: + 7V)]i<,-,<Ar = 

J2 det[4^(M. - z + iV, Aj - J + A^)]i<.j<jvdet[/j'^'°)(/i, - t + iV)]i<.j<Ar, 

which implies the proposition. 

For a — ±1/2, define the matrix T"^ ^ on J^r x Jtv by 

T- M A) = det[/r(A.. - ^ + ^, A, - , + ^^^,,,^^^^^^^^±^1^. 

<i™SO(2Af+l/2+a) M 



□ 



Proposition 3.1 suggests a Markov Chain with state space lN,a with transition 
probabilities given by ^ . 

Proposition 3.2. For any ip G C^[— 1,1], the rows of T"^ ^ sum to (p(l)^. In 
particular, if ip{\) — 1, then the rows of T'^ ^ sum to 1. 

Proof. Let zi, . . . , z^r be complex numbers on the uni t cir cle and set Q = {zi 
z^^)/2. Using the notation and statement of Lemma 



2.3 



-1/2) 



l<i,j<N 



det[g{fii -i + N,Xj - j + N)] i<t,j<N 



det 



,(a,-l/2) 



(0) 



l<i,j<N 



By (111, for any k,l>0, 



Therefore 



(q -1/2) 
\j-j+N 



(0) 



l<ij-<7V 



det [I^{fi,^i + N, - .J + N)]^^^^ 



det 



(1,-1/2) 



l<i,j<N 
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or equivalently, by Q and (|5|, 



^SO(2N+l/2+a) {zi,...,ZN)det [I^ {n^ - i + N, \j - j + N)]i<,j<N 

= XsO(2W+l/2+a)(^l' ■ • ■ ' Zn)'P{Ci) ■ ■ ■ ViCN)- 

Taking zi, . . . , zjy = 1 shows that the rows of ^ sum to Lp{l)'^ . □ 

Proposition 3.3. If ip{x) ^ po+pix with po > pi>0, then each entry ofT^^ 
is nonnegative. The same holds if ip{x) = e*'-^^"'^-', where t > 0. Additionally, 
the diagonal entries of T^'^^^j'^ are bounded below by 

( R\{R+ -pi/2)-RLiR- -pi/2)\ 



(17) 



where 

R± = 



Po ± Vpo^Pi 



Proof. First consider the situation when ip{x) = po+pix. Note that by ([8|-( 10 1, 
I^ikJ)^Oif\k~l\>l. 

If /ii < Ai — 1 for some i then < A; — 1 for k > i and I < i, which implies 
that I^{fj,k,)^i) = for such k,l, and thus the determinant in question is 0. If 
fii > Xi + 1 then /i/j > A; + 1 for k < i and I > i, which means laifJ'k, A;) = 0, 
and the determinant is again. So det[/^] is zero if |Ai — /Uij > 1 for some i. 
Hence, it remains to consider the case when |Ai — /J,i| < 1 for all 1 < « < A^. 

Split + iVjili into blocks of neighbouring integers wth distance be- 

tween blocks being at least 2. Then it is easy to see that det[Ia{fJ.i — i + N, Xj — 
j + N)] splits into the product of determinants corresponding to blocks. It 
suffices to show that the determinant corresponding to each block is nonneg- 
ative, so assume without loss of generality that — i + -/Vjfli is one such 
block. In other words, assume that all /ij are equal. Then there exist m and n, 
1 < m < n < iV such that /i^ = A; — 1 for I < i < m, and /i.; = Xi for m < i < n, 
and = Ai -|- 1 for n < i < N . The determinant is the product of determinants 
of three matrices. We shall examine each of these matrices. 

Because I^{k,k ± 1) > 0, the matrix parametrized by 1 < i,j < to is 
triangular with nonnegative diagonal entries, so has nonnegative determinant. 
Similarly, the matrix parametrized hy n < i,j < N also is triangular with 
nonnegative diagonal entries. It remains to consider the matrix parametrized 
by TO < I, j < n. 

For now, assume A„_i =/= 0. Then this matrix is tridiagonal, with po in the 
diagonal entries and pi/2 in the subdiagonal and superdiagonal entries. If this 
matrix has size r x r, let Dr denote its determinant. Then Dr satisfies the 
recurrence relation 

9 2 
P p 

Dr = PoDr-1 - ^Dr^2, Di = po, 1^2 = Po " ^ • 
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Solving this explicitly yields 

r+1 / -\ r+l^ 



Dr = 



Po 



VpI -pI] (Po-y^ 



VpI - pI 



which shows that Dr is nonnegative. 

If A„_i ~ and a = 1/2, then the entry in the rth row and rth column is 
Po ~Pi/2 instead of po- The other entries are pQ,pi/2, and 0, as before. In this 
case, the determinant is 



,2 

'^r-2 




Po - Vpq-Pi 



VpTpI^ 2 



Note that Ar is positive on the interval (pi/2, oo), which contains {po+\/Po ^Pi)/2, 
and is negative on the interval (0,pi/2), which contains (po — \/Po ~ P'l)/'^- 
Therefore the above expression is nonnegative. 

If A„_i = and a = —1/2, then the only modified entry is in the rth row 
and (?■ — l)st column. It equals pi instead of pi/2. In this case, the determinant 
is 

2 2 
Pi_ P\ 
PoDr-1 - -^Dr-2 = Dr - ^A_2- 

We have already shown that Dr > {pi/2)Dr-i, so therefore Dr > (Pi/4)-Dr-2- 
So the above expression is also nonnegative. 

For the last claim in the proposition, note that Dr > Dr — (j3i/2)^_Dr-2 > 

Dr-{pi/2)Dr-l = ^. 

Now let ip{x) = e*(^-i). For n > 2t, let (/3„(a;) = (l + ■ By Lemma 



3.4 



below, T^;^ = T^^'a""^^''' ■ ■ ■ T'wta*''"^^^"- We just showed that all entries of 
each are nonnegative. Therefore all entries of T^"^ are nonnegative. 



By Lemma 3.5 ^{fj,,\) equals lim„_>oo ^'^'^(Mi •^); so all entries of ^ are 

□ 



nonnegative. 

Lemma 3.4. //v?i,¥'2 e C^-l, 1], then T^J^ = T'^n^T'n^- 

Proof. This is a straightforward computation using the Lemmas 2.1 and 2.2 □ 

Lemma 3.5. Suppose {(pn)'^^=i,^ € C^[—l, 1] and (fn converges to (p uniformly 
on[— 1,1]. Then for any \, p £ Jin , T'^\{^,\) converges to T"^ ^{pL,\). 

Proof. Since (/?„ converges uniformly to Lp and Lp is bounded, the dominated 
convergence theorem implies that each converges to Since 

T^"^(/i, A) is continuous in the variables I^"'{pi — i + N, Xj — j + N), it must 
converge to ^{p, A). □ 
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3.2 Generalities on Multivariate Markov Chains 



Recall that in section 2.2 we explained that the measures P"^ ^ generalize to P"^. 
We just constructed a Markov Chain that maps P'^ ^ to P'^ ^, so it is natural 
to expect a Markov Chain that maps P'^ to P'^. Our goal now is to extend 

-I- to stochastic matrices on finite paths that map P'^ to (or rather their 
projections on finite paths). This section describes a general construction from 
U which builds a Markov chain from smaller ones. In the next section, this 
construction will be applied to our case. The original idea for bivariate Markov 
chains goes back to f^D]. 

Let Si, . . . ,Sn be discrete sets. For 1 < fc < n, let be a stochastic matrix 
with rows and columns indexed by Sk- For 2 < fc < n, let A^_j^ be a stochastic 
matrix with rows indexed by Sk and columns indexed by Sk-i- Assume these 
matrices commute: 

Ati := AtiTfc_i = r,Ati. 
The state space for the multivariate Markov Chain is 

oo 

4"^ = {(a;i,...,x„) e5i X ... X 5„ : Y[A''k_^{xk,Xk-i)^0}. 

k=2 

Write Xn = {xi, . . . , a;„), y„ = {yi, . . . , ?;„) G 5^""*. The probability of a transi- 
tion from Xn to Yn is 

.-fr Tk{xk,yk)At-iiyk,yk-i) TT^fe . 
k=2 '-^k~i\-^k, yt-i) j,^2 
0, otherwise. 

Let T denote this matrix of transition probabilities. One could think of T 
as follows. 

Starting from X — {xi, . . . ,x„), first choose t/i according to the transition 
matrix Ti{xi,yi), then choose y2 using ^ which is the condi- 

tional distribution of the middle point in the successive application of T2 and 
A^ provided that we start at X2 and finish at yi. Then choose y^ using the 
conditional distribution of the middle point in the successive application of T3 
and A2 provided that we start at X3 and finish at y2, and so on. Thus, one could 
say that Y is obtained by the sequential update 

The next proposition will be used later. 

Proposition 3.6. Let mn{xn) be a probability measure on 5„. Consider the 
evolution of the measure m„(a:„)A5^_]^(a;„, . . . Af{x2, xi) on iS^"-* under the 

Markov chain T, and denote by {xi{j), . . . , Xn{j)) the result after j = 0, 1, 2, . . . 
steps. Then for any fci > ^2 > . . . > fc„ > 0, the joint distribution of 

{•^n (1) ; • ■ ■ 1 •^ni^n ) ; *^n — 1 (^n) i *^n— 1 (^n ~t~ 1 ) • • ■ ; -^n— 1 (^n — 1 ) : 

a;„-2(fcn-l), • • ■ ,X2{k2),Xi{k2), . . . ,Xi(fci)) 
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coincides with the stochastic evolution of nin under transition matrices 

kn /C„_i — fc„ fcl — fe2 

Proof. See Proposition 2.5 of [8]. □ 

Let L be the linear subspace of spanned by elements of the form 

mn{xn)-An-ii^n,Xn-i) ■ ■ . A^(a;2,xi), where to„ is a summable function on 5„. 
Then T can be thought of as a bounded linear operator of L. Similarly, T„ is a 
bounded linear operator on 

Lemma 3.7. With the notation from above, 

\\T-id\\L = ||r„-M|hi(5„). 

Proof. By definition, 

IIT Ml IIt/-/||l 

\\l-ld\\L — sup r— , 

feL II/IIl 

||T„ - = sup r r . 

m„e/i(5„) ll'™n||;i(5„) 

If / G L, then / must have the form rn„(a::„)AJJ_]^(a;„, a;„_i) . . . Ai{x2, Xi) 
for some m„ G Since the matrices A^^"'^ are all stochastic, ||/||l = 



ll"in||ii(5„)- By Proposition|3 6l T/ = (r„m„)(x„)A;^_i(a;„, a;„_i) . . . Af (a;2, Xi), 
which imphes ||T/ — = ||'J„m„ — m„||;i(5^). Thus, the lemma holds. □ 

3.3 Markov Chain on Multiple Levels 

We need an implementation of the stochastic matrices A^ -^. Define the matrix 

T^'^ on Jat + X Jat^- by 

A) = 5"^""^^"^ ' det[r(^, - . + iV, A, - , + N)]^^, 
dimso(2Af+i) M 

where 

T{x,y)= <^2, x>2;>0, . 
[o, x<y. 

Also define the matrix T^Ci + on In.- x In-i,+ by 



A) = ' det[<^(M. -^ + N,X,-J + N-r 



where 



1, x> y, 
0, a; < y. 

Recall the definition of x in Section 12.21 
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Lemma 3.8. For any (A,^) G JAr-i,+ x In-, 
det[0(M. -i + N,Xj~j + N ~ 



For any {X,fJ,) G In,- x In, 



Proof. The argument is standard, see e.g. Proposition 3.4 of [II . The proof of 
the second formula is exactly the same. 

□ 



Proposition 3.9. The matrices T^'^ and T^'' 



are stochastic. 



Proof. First let us show that T^'T is stochastic. Taking dimensions of both 
sides of (12) yields 



Therefore 



By Lemma 3.8 T^'^ has nonnegative entries, so it is stochastic. 



Now we will show that T^'_^^ is stochastic. Taking dimensions in (13 1 
yields 

dimso(2A.) ^^ - E dim50(2^-i) A. 

Therefore 

rw-"i,+ (M,A) = 1. 

A(EJjv-1,1/2 

The nonnegativity also follows from Lemma |3.8| □ 



The matrix elements of Tj^'J*" and T, 



'^N-l,+ 



are cotransition probabilities of 
the branching graph J, see e.g. [32]. Recall that in Section 3.1 we defined 
stochastic matrices Tn. ±1/2- We use the notation Tn.± for convenience. 

Proposition 3.10. Assume (/?(1) = 1. For any N > 1, we have the following 
commutation relations: 



Ttip rpN,-i- rpN , + rp^ 



^N,-^N-1,+ — ^N-1,+^N-1,+ - 
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Proof. We start by proving the first relation. By Lemma 2.1 



rr.N.+ l ^^ dim A 

X ^ det - i + N,v.j - i + N) 

oo 



dim A 
dim/x 



det 



^ -i + N, x)T {x, Xj -j + N) 



.x=0 



l<iJ<N 



Similarly, 



dim A 
dim II 



det 



J l<iJ<N 



_x=<3 

We thus need to check that 

C30 

^ T(^, - z + TV, x)/!,/2(.T, A, - j + iV) 

oo 

= ^ /r/2(A^* - I + A^, x)r{x, Xj - J + N). 



x=Q 



x=0 



By applying Lemma 2.7 to the right hand side and Lemma [2^ a) to the left 
hand side, one sees that both sides are equal to 



Now we prove the second relation. Expanding det[0]^ along the TVth column 
and using Lemma |2.1[ we obtain 

(18) 



,+ 

dim A 

dim /i 



N 



II(-l)'^"^' det -i + N,v,- ] + N ^ l)]i<.^fc<iv 



X det[/^/2(i^^ -i + N ~\,Xj - j + N - \)]i<ij<N-i 



l<j<jV-l 



(19) 
(20) 



dim A 

dim/i 



N 



^(-l)^-'=det 



fe=i 



^ (/.(^, -i + N, x)I'C^^{x, Xj-j + N- 1) 



x=0 



l<i^k<N 
l<j<N-l 



dim A 

dim li 



det 



4>{^^^ - * + iV, r)I^,2{r, A, - J + TV - 1) 



.r=0 



(21) 
(22) 



J l<i,j<N 
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where it is agreed that all matrix elements in the A^th column (j = N) of 22 
are equal to 1. Similarly, 



dim A 



dim fi 



det 



E ^-i/2(a^^ -i + N, r)0(r, A, - j + TV - 1) 



r=0 



(23) 



l<ij<JV 



2.5 », the Nth column of J23]) equals ¥3(1)J^^1{^'^^/^^(1) = 1. 



By Lemma 

Therefore the iVth columns of (22 1 and (|23p are equal 



By Lemma 2.4 'b), for j ^ N, the (z, j)-entry of (22 1 equals 

''•-'/'\{xMx){l - + xy'/'dx, (24) 



1 , (-1/2,-1/2). X 

^ -'fii-i+N W_ .(1/2,- 



1-a; 



'A, 



By Lemma 2.5 'b), for j < N, the (i,j')-entry of (23 1 equals 

»i 



1 _ |(-l/2 -1/2) 



{x)v{x))}^^ll^r^%%{x){l x)-^'\l + xy/'dx. (25) 



Their difference only depends on j, so for j ^ N, 



jth column of (23 1 = jth column of (22 1 + [(25 1 - (24l](7Vth column), 



so the matrix in ( 23 1 is obtained from the matrix in ( 22 1 by elementary column 



operations. This means that their determinants are equal. 



□ 



Following the notation of section |2.2| let us denote the set of finite paths in 
of length (27V - 1/2 + a) by 



j(W),a |;^(l),-l/2 ^ ;^(l),l/2 



^ . . . ^ 



Using the construction in section 



3.2 



the stochastic matrices _|_ , T^'^ , T^^^ ^ 
allow us to construct a Markov Chain on J^^^'". Define the matrices A]^'^ and 

^N-i,+ by 

AA'',+ rpN, + rp^ rp^ rpN ,+ 



Define the matrix A'^^,^ with rows and columns indexed by elements of 
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j(JV),l/2 by 

<.i/2(^^'''"'^' ^ . . . M^'^^"'/' ^---^ 

C\(l),l/2 (l),l/2^yl,+ C„(l),l/2 (1) -1/2N 

^ A};^(A(i)4/2, ^(1)^-1/2) 



. . X 



A;^'+(AW,i/2^^(JV),-i/2) 



if A^;1(aW'1/2,^W.-i/2), A^+i'-CAe^+i^-i/^^/^e^)'!/^) ^0, 1 < < TV 

0, otherwise. 

Simiarly define with rows and columns indexed by elements of J^^'''^^/^ 

by 

(XW■^/^ (i),i/2 (1) -1/2N 

''-^ ^ A};l(A(i)-i/2,^(i),-i/2) 

r;^__(A(^)-i/2^ /iW-i/2)r^-^^(^(^).-i/2^ ^(^-i).i/2) 

■■■ A^Li,+(AW.-V2,^(A^-i)a/2) 

if A^;+(AW'l/2,/xW-l/2)^^fe+l,-(;^(fc+l),-l/2^^W,l/2) q, 1 < < TV 

0, otherwise. 



By the construction in section 3.2 ±1/2 stochastic for tp{x) = l—pi+pix, 

< Pi < 1/2, and ip{x) = e*(^^,t > 0. 

For any E{x) £ - lA] such that £^(1) ^ 0, let pC^)'-^ be the (possibly 
signed) measure 

Pjv,a(AW^'^) . . . Tf;-(A(2)-i/2^ A(i).i/2)ri.+(A(i).i/2^ A(i)-i/2) 



3.1 



on 

j(w),a^ where P^v, 

a is as in Theorem 2.8 Proposition 3.6 and Proposition 

imply the following. 

Proposition 3.11. Let E{x) e C^[-l,l] such that E{1) ^ 0. Consider the 
measure PW-'^ onjW'". After one step of the Markov chain A^ ^, the resulting 
measure on IS 

Piv,a(AW^'')...Ti2;+(A(2)-V2^A(i)^i/2)Tii;+(A(i)^i/2,A(i)'-i/2), 
where Pnm is defined from the function E(x) = Lp{x)E{x). 



3.4 A Continuous-time Markov Chain on Multiple Levels 

Define a matrix Qn a on J'^^-''" x J^^)-'' as follows. Let us explicitly write 

QiV,a(A«'-l/2 ^ A(l).l/2 ^ ... ^ AW,a^^(l),-l/2 ^ ^(l),l/2 ^ ... ^ ^W,a)^ 

There are three cases to consider: 
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Case 1. This occurs when there exist (no, uq) < {ni,ai) and k < no such that 
the numbers — 1, A^" are all equal for (no,ao) 1^ {n*,a*) < (ni,ai). 

Furthermore, /il"-*'" = other h,a,l. 

There are two subcases: 

Case la. When case 1 is satisfied and ao — —1/2 and X^^°^'°'° = Q. 

Case lb. When case 1 is satisfied and case la is not satisfied. 

Case 2. This occurs when there exist (no,ao) ^ (ni,ai) and k < uq 
such that the numbers a!" If, , If , . ^ + 1 are all equal for 

(no, flo) < (f^*, a*) ^ ('T-i, oi)- Recall that (i(ni, ai; n-o, ao) = |2ni+ai — 2rio— ao|. 
Furthermore, — for all other h,a,l. 

Case 3. This occurs when the two paths A and fi are not equal and neither 
case 1 nor case 2 is satisfied. 

When case lb or case 2 occurs, the corresponding element of QN,a is 1/2. 
When case la occurs, the corresponding element is 1. When case 3 occurs, the 
corresponding element is 0. The diagonal entries are defined so that the rows 
of QN,a sum to 0. 

Under the map £23, the cases can be described more easily. Let {j/™} = 
/!<g(A) and {2;™} = C<x){fj,). Case 1 occurs when there exist mg < mi and 
k < [^] such that 

- 2 = = »+i - 1 - 2 = yp+^ - 1 = . . . 

= z["i - (,ni - ,no) - 2 = - (mi - mo). 

Furthermore, z™ = y[" for all other l,m. 

Case la occurs when case 1 is satisfied and y™° = 0. 

Case lb occurs when case 1 is satisfied and case la is not satisfied. 

Case 2 occurs when there exist toq < "^1 and k < [^] such that 

z™o + 2 = = z[-+i + 1 + 2 = + 1 = . . . 

= ^r+mi-mo + ("^1 -mo) +2 = + ~ "^o)' 

Furthermore, z"' = for all other l,m. 

Case 3 occurs when {y™} ^ {-z™} and neither case 1 nor case 2 is satisfied. 

It is not hard to see that QN,a is the generator of the continuous-time Markov 
Chain defined in Section [T] In general, if Q is a matrix with countably many 
rows and columns such that its rows add up to 0, its ofi^-diagonal entries are non- 
negative, and its diagonal entries are uniformly bounded, then there is a unique 
continuous-time Markov chain with Q as its generator (see e.g. Proposition 2.10 
of [T]). In words, this Markov chain satisfies 

• In state i, a jump takes place after exponential waiting time with param- 
eter Qii- 

• The system makes a jump to state j with probability ^Qij/Qa- 
We aim for the following: 
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Theorem 3.12. Let tpt{x) = e**^^ ^\ Let P^^^'°- he the (possibly signed) central 
measure on J'-"'^-'''' corresponding to some E{x) e C"'^[~l, 1] satisfying E{\) ^ 0. 
Then 6*0"." • = • p(^)^\ 

Proof. This theorem reUes on the following proposition. It can be found as 
Theorem 9.6.1 in [24]. 

Proposition 3.13. Let {A(t) : t > 0} &e hounded linear operators on a Banach 
space B such that A(s + t) = A(s)A(t) for all s,t > 0. Lf lim \\A(t) - /|| = 0, 

then there exists a bounded linear operator Q on B such that A{t) = e**^ for 
t > 0. 

Each Al'^ ^ is a linear operator on /^(J'^^'-'^). Since its matrix is stochastic, 
it is a bounded operator. 

Let L be the linear subspace of ?"'^(J^^^'°) spanned by all measures corre- 
sponding to functions F G C^[— 1,1] satisfying f (1) ^ 0. Define the Banach 

AVt AVs _ AVt+s T 



3.11 



space B as the completion of L. By Proposition 
By continuity, the same holds on B. So it suffices to show that — /|| ^ 

By Lemma 3.7 it equivalent to show that HT^'^j — /|| ^ 



More precisely, it must be shown that 



lim sup \T'^\a{\l^)-5xA^Q. 



Aej(«) 



Since T"^^^ is stochastic, it is equivalent to show that 

lim sup (2-2T;^'JA,A)) = 0. 

So it sufllces to show that lim inf T)^* (A^^)'", A^^^'") = 1. 

We prove that T;^4(AM'±, A('')'±) > exp(-i(r + 1/2)). Let (/?|"' = (l + t(x- 
l)/n)". By Lemmas |3.5| |3.4[ and Proposition |3.3| respectively. 



r,^i(AM^±,A('-)'±) = lim r;^| (AW^±,AM'±) 

' n— *oo ' 

= lim (r;+*("~i'/")"(AM'±,AM'±) 

n — *oo 

> lim [T,^+*("-^)/"(A«'±,AW^±)j" 

From Proposition |3.3[ 

yl+t(a;-l)/n^^(r),± ;^(r),±^ > ( ^-Pl/^) " RL{R- ~Pl/2) 

where 



t _ t _Po± Vpq-Pi 



Po = 1 , Pi = -, R± 

n n 
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Since i?_ < Pi/2, 



T,^i(AM>±,AM'±) > lim 

' n — >oo 

Finally, notice that as n — > oo, 



R7>{\/pi-pi] =(1--) 



-rt 



We have just shown that •P(^)'" = e*'Q.p{N).a g^j^^g gj^jgi^ ^j^g 

proof, we show that Q — ^A^*^|t=o — QN,a- Since we only need to calculate 
Q up to terms of order O(t^), we can replace '^t{x) = e*'^^^' with \ — t + tx. 

The problem now is to calculate A^'^ up to terms of order 0{t^). There are 
three cases to consider: when all the particles on the mth level stay still, when 
one of the particles on the mth level is pushed by a particle on a lower level, 
and when one of the particles on the mth level moves by itself. As an example, 
consider particles on the (m, 1/2) level when one of them is pushed. 

The expression that needs to be calculated is 



A:::;i/,%(A('")'i/2,^M,-i/2) ■ 



(26) 



Assume that ^ Since /i('")'-i/2 ^ ^(™).i/2^ this implies 

that A^'"''!/^ ^ ^(")'i/2. Similarly, A^™) "^/^ ^ ^(m),-i/2^ ^Yach means one of 
the particles on the (m, —1/2) level is pushing a particle on the (m, 1/2) level. 
Conversely, if the fcth particle on the (m, — 1/2) level is pushing a particle on 
the (m,l/2) level, then ^i'^^'>'-^'^ > A^^"^''/', so ^M-i/2 ^ aM'1/2. 
The transition probability on the (m, 1/2) level is (because of (Is])) 



^ /\(™).l/2 „M,l/2^ t dimgo(2„,+l) , ^^^2^ 

Jm,l/2(A' ^' ' ' j = + U[t ). 

2 dim5o(2m+i) A("'^V2 

Furthermore, 

. m,l/2 / ,(m),l/2 „(m),-l/2N „ dimgo(2^) ^('") "1/2 

"-^/^^^ ^"dim50(2™+i)A*(")-i/^' 

A:;i/,%(A('")'i/^M^'")'-V2) 

= C:i/^^/2(A(™)' A(™)-i/2) • ^,„,_l/2(A(")-l/^ A*^'")-^/^) + o(t2) 

_„...W,-i/2 /"^50(2„.)A('")'-V^ ^ M/(^(r)-'^/^)dimgo(2„)^('").-V^ 

^dimso(2™+i)A(™).i/2 2 vt/(A(^')'-i/2) dimso(2„.) AM.-1/2 + 
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Therefore equation (26 1 equals 1 + 0{t). 

Similarly, when all the particles on a level stay still, the contribution is 
1 + 0{t). When a particle against the wall moves, the contribution is t + O(t^). 
When a particle not against the wall moves without being pushed by a particle 
on a lower level, the contribution is t/2 + 0{t^). □ 



4 The Correlation Kernel 

Theorem 4.1. For any w € f2 with parameter Pi < I, the point process 

is determinantal. Denote its correlation kernel by _ftr"(ni, oi, Si; n2, a2, S2). // 

(ni, ai) > (^2, 02), then K'^ {ni, ai, si; 712, 02, S2) equals 



1 



X 



{x - (1 - x)''^ (1 + x)-^/^dudx 



(u — X — u 

// (ni, ai) <l (^2, 02), then K'^{ni,ai,S2;n2,a2,S2) equals 

{x - (1 - x)"^ (1 + x)-^/^dudx 

X "~ r . 

[U — X — u 

The u-contour C is a positively oriented simple loop that encircles the interval 
1, 1] but does not encircle any zeroes of E'^ . Recall that the functions -E'^, Js 
and W were defined in Section\2.1\ 



Remark. The case /3i — I can be obtained by the limiting transition /?i ^ 1 
from Theorem 14.11 



Corollary 4.2. With the definition of the particle-hole involution A given in 
Appendix \A\ i/ (^2, 02) oi), then K'^{ni,ai,si; 712, a2,S2) equals 

{x - l)"i (1 - x)-^! (1 + x)-^/^dudx 

{u — 1)"2 X — u 
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If{ni,ai) < (^2,02), then K'^{ni,ai,S2;n2,a2,S2) equals 



{x ~ (1 - x)"^ (1 + x)-^/^dudx 



{u - x - u 

Proof. This result follows from the orthogonality relations 



dx 



for a = ±1/2. Note that in Theorem 4.1 the two cases are (^2,02) < (ni,ai) 



and (ni, ai) < (712, 02). Here, the two cases are (7^2, 02) <i (rii, ai) and (rii, ai) < 
(712,02). □ 

This proof uses Theorem 4.2 from [7 , which we will describe in the next 
subsection. We will alter the notation to make it more convenient later. 



4.1 Determinantal structure of the correlation functions 



Theorem 2.8 and Lemma |3.8| imply that the measure of a finite path A 



N 

= const • 

ra=l 

X det 



det 



("),i/2 , , „ x(")-i/2 J. , „^ 



l</c,i<^ 



X det 

where 6 and T were defined in Section 



i<fc,;<7v 



(27) 



3.3 



and J^'"" was defined in the statement 
of Theorem 2.8 If a = —1/2, then the final determinant with the T does not 
occur. 

Recall that we have set aL" ^•''^^^ to be equal to zero, so An ^ ^ ri -\- 
n ~ 1 = —1. We will refer to —1 as a "virtual variable," or "virt." 
Set 

TJ/(a.-l/2)/„\ l-l 

K-ii-^)^- — / E^{x)^t-'^^Hx)ix-^)''-\l~xn^ + x)-'/^dx. 

^ J~i 

(28) 

Observe that Span{/i, . . . , /jv} — Spanl^fi, . . . , ^'tv}- Thus, if we replace /; in 
(27) with "^N-i, the measure is not going to change. 
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Let * denote convolution. More explicitly, 

(/ * 9)ix, y) = Y^ fix, z)g{z, y), (/ * h){x) ^ ^ f{z)h{z, x). 



z>0 



For (ni, fli) <i (n2, 02), set 



i(ni,ai),(n2,a2) 



ry 



"2 — "1 — 1 



if ai = -1/2, 02 = 1/2, 
if ai = 1/2, a2 = 1/2, 
if ai = 1/2, a2 = -1/2, 
if ai = -1/2, 02 -1/2. 



For (ni,Oi) \> (n2,a2), set 0("i,ai),("2.a2) ^ g. 

Let M be the N x N matrix with entries 



kl 



- N- 

"^N-i * * * (t))^-''{virt), if a = -1/2. 



(virt), if a = 1/2, 



For k<N, define = ^-^'^^^ * 0(n,ai),(jv,a)^ 



Theorem 4.2 from j7] also says that if M is upper triangular and invertible, 
then there exist functions ^^^'^l{s) such that 



if 02 = 1/2 
if a2 = -1/2 



• {*J'n-fc('S)}fc=i....,n is a basis of the linear span of 

\{((0*r)"-'=*(/.)(s,mrt)},=i,...,„, 

• For < i, j < n — 1, 

s>0 

The formula for the correlation kernel is given by 

"2 

K{nuai,sr;n2,a2; S2) = -</>("-'^^)'("-"^)(s2, si) + ^ vl/«;:_-^(si)<i>::^L"^(s2). 

fe=i 

(29) 

In section 4.2 we pro ve th at M is upper triangular. In section |4.3| we 



calculate ^' }. In section 

ni — ft! 



4.5 



4.4 



we calculate ^ *^'':°fc(si)*;^2:!fc(s2). In section 

we add all these 



k = l 



we calculate (^*^"^'"i^'("^'°^^(s2, si). Finally in section 
expressions together to get the expression in Theorem |4.1 



4.6 



33 



4.2 The Matrix M 

Lemma 4.3. The matrix M is upper triangular and invertible. 
Proof. By definition, 

'Y,^^-^i{s){T * cf>)^-''+\s, vtrt), a = 1/2, 

s>0 

J2 K'M ^ct^HT* 4>)''-\s, virt), a = -1/2. 



kl 



Define 



{T *(f))'^-''+^{s,virt), a = 1/2, 
(l)*{T *<j>)'^-''{s,virt), a = -1/2, 



so that Alki = (^'^'°;, g^"^^). Tiie definitions of T and (f) imply that 5^°^(s) is a 
polynomial of degree 2N — 2fc + 1/2 + a. 
Define 

where the integration contour is a positively oriented simple loop around u = 1 
that does not contain any zeroes of E'^{u). (If /3i < 1, then E'^{u) has no zeroes 
in [—1, 1], so the contour always exists). The integrand has a pole only at u = 1, 
hence 

,N.i 



1 d^'^ 1 



Clearly, $jY-fe('^) ^ polynomial of degree 2N — 2fc + 1/2 + a. 

Since g^"^ and are polynomials of degree 2N — 2fc + 1/2 + a, there 

exists an invertible upper triangular matrix A such that 

N 

9k - 2^ Akm^N-m- 
m—1 

Therefore 

JV AT 
m—l m—1 

Finally, by Lemma [2. 2 1 

SO Af = j4, which is upper triangular and invertible. □ 
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4.3 Calculating 

For ease of notation, let E denote E'^ in the remaining sections. 
The purpose of this section is to prove the following: 

Theorem 4.4. For I < ni, 



|4A(ai,-l/2)(g) 



E{x)^i''''-^/^Hx){x-iy''-\l-x)'''{l+x)-^/^dx. 



(31) 



For I > ni, 



|y(ai,-l/2)(g) 



1 E{x) - E{1) - E' {l){x - I) 



'(1) 



(/-m-l)! 



{x- 



[x - 

X J('^i^-i/2)(x)(l - xY'il + x)-^'^dx. (32) 



Proof. Start with the proof of ( |31[ ) . It will be done by induction on d(ni ,ai;N, a) 
2{N - m) + a ~ ai. When d{ni,ai; N,a) = 0, then ([3l| is true by Now 
assume that (31 1 holds whenever d{ni,ai; N,a) = m. Either ai = 1/2 or 
ai 

ing Lemma 



-1/2. If ai — 1/2, then convoluting both sides of (31 1 by T and apply- 
with T{x) = E{x){x - l)"i^' shows that pl| holds for ni and 
ai = —1 /2. If oi = —1/2, then convoluting both sides of (|31[) by and applying 



Lemma 



2.5 



1 



1. 



» with T{x) = E{x){x - 1)"!-' shows that (jSlj holds for ni 
and ai = 1/2. Either way, (311 must hold whenever d{ni,ai; N,a) = ni 

Now on to the proof of ( 32 ) . It also will be done by induction on (i(ni , ai ; iV, a) . 
The base case occurs when I — ni — 1, and ai — 1/2. We just proved that (31 1 
holds when ni = I and ai = —1/2. Convolute both sides of (31 1 by and apply 
Lemma 2.5 ^b) with T{x) = E{x). This proves the base case. 

Now assume that (32) holds for some ni and ai = —1/2. Convolute both 



sides of ( 32 ) by 



Apply Lemma 2.5 'b) by setting 



T{x)^ 



E{x) - E{1) - £;'(l)(a;- 1) 



(a; - l)'-"i 



(i-m-i)! 



(x-1) 



This shows that ( 32 1 holds for ni — 1 and oi = 1/2. Now assume that ( 32 ) holds 



for some ni and Oi = 1/2. Convolute both sides of (32 1 by T. By Lemma 2.7 
with 



T{x) 



E{x)~E{l)-E'{l){x-l) 



(i-m-i)! ^ ' 



1 



{x - l)'-"i 



(32) also holds for rii and ai = —1/2. 



□ 
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4.4 Calculating 

Define for 1 < fc < n < A^, 

1 r 1 s du 



2TTi J E{u) " ' '(u- 

where the contour contains the interval [—1,1] and does not contain any zeroes 
of E{u). Note that this agrees with (30 1. 

a , 
-k\ 



Lemma 4.5. (a) {<i>^!°j,(s)}fc=i^....„ is a basis of the linear span of 



I {(T * (0 * T)^-" * ^)is, virt)}k=i,...,n, ^fa^ 1/2 
* r)"-'=) * 0(s, mrt)}fc=i,...,„, z/ a = -1/2 

/6/ For < i,j < n - 1, 

Proof, (a) *i>^i!^i,(s) only has a pole at u = 1, so it equals 

1 ji^'^^/^'Hi 



(n - fc)! da;"-*^ E{u) '"=!' 

which is a polynomial in s of degree 2n — 2fc + 1/2 + a. Also, (T * (0 * j-yi-k ^ 
0)(s, wirt) is a polynomial in s of degree 2n — 2k + 1 and {{(p*Ty''^'^ *4>){s, virt) 
is a polynomial in s of degree 2n — 2k. This proves (a), 
(b) By Lemma [2. 2[ (the contour below contains [—1,1]) 

du 



s>0 ^ ' 



□ 



Proposition 4.6. In the expressions below, the u-contour is a positively ori- 
ented simple loop that encircles the interval [—1, 1] but does not encircle any 
zeroes of E. 

For any ni > n2 > I and Si, S2 G Z>o, we have 



k=l 



(x - l)"i (1 - a;)"i(l + x)-^/^dudx 
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If ^ l£ n-1 < n-2, then 



fc=i 



H^(ai,-l/2)(^^) rl j(ai,-l/2)j(a„-l/2) 



TT 



-1 
X 



Eix)-E{l) 



^(Ti2-ni-l)^2) 

(n2 - ni - 1)! 



M/(--v2)(,^) 1 /I / -E{i) ... - 1)---^-^ 



2TTi 



x4^--l/2)(^)j(a.,-l/2)(^)(^_l)r 



(1 - x)''^{l + x)-^/'^dudx 



X — u 



+ 



py-(ai,^l/2)(g^) 1 



E{x) 



j(a.-l/2)(^)j(a.-l/2)(^) 



(x - (1 - x)''! (1 + xy^/'^dudx 

K r . 

[u — X ^ u 



Proof. First assume that ni > n2. Then the left hand side of (33) equals 

^(ai.-l/2)(^^) 1 



TT 27ri 
The expression 

^(ai,-l/2)(g^) 1 
TT 27ri 



(a; - 



u - 1 

X- 1 



(l-.T)''i(l + a;)"i/2^x 



j(a.-l/2)(^)j(a.-l/2)(^) 



E{u) ' ' 

(a;-l)"i {l-xy^{l + x)'^^^dx 



X — u 



(u- V^;- 1, 
has residues only at u — x, so equals 

^ ^(ai.-l/2)(^^) j(^i-l/2)(2.)j(a2,-l/2)(^)(^_l)ni-n.(^_^)ai(^^^)-l/2^^^ 

(34) 

Hence (33 1 holds. 

"1 

Now assume m < n^. ThenJ^Kl-li-'^Xl-li''^^) equals 



^(ai.-l/2)(^^) 1 



Eiu) 



j(--V2)(,)j(a.-l/2)(^) 



1 - 



u- 1 

X- 1 



x ~ u 



(35) 
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To evaluate ^ Kl-ki-'^i)Kl-kis2), we first evaluate 

k—'rii-\-l 



E{x) - E{1) E'{l){x - 1) - ... - - 1)''"^'' 



E 

fc=ni + 1 



It can be rearranged as 



Each sum is a geometric series, which can be explicitly evaluated. After simpli- 
fying, we get 

(y _ \ \ni-n2 ( 



X — u 



(n2 — ni — 1)! x~u \ x—1 



Therefore ^ ^:iL°Usi)<i>:^"^(s2) equals 

fc— ni + 1 



i^(--v2)(,^) 1 .1 . - E{i) - ... - - 1) 



TT 27ri J^iJ E{u) 

X j(^--l/2)(:,)j(-2^-l/2)(^)(^ _ (1 - x)"^ (1 + x)" ^^do; 

X — u 

^(a„-i/2)(,^) 1 .1 . E{x) - E{1) - ... - ^g^^gW (u - 



TT 27ri J Eiu) 

X j(''--V2)(a;)j(--'V2)(y)(^_ l)™l-n.(l^^)!l(l±^)l!^. (36) 



The first term in (36 1 has residues only at m = x, so simplies to 

H^(ai.-l/2)(^^) .1 j(ai,-l/2)j(a„-l/2) 



-1 

- - {x-iy''-''^{i-xY'{i+x)-^/^dx. 

(n2-ni-l)! y 
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Adding ( 35 ) and ( 36 1 and rearranging finishes the proof. □ 

4.5 Calculating ^("i.^O.C'^Z'^a) 
For (ni, fli) < (n2, a2), set 

r(rii,ai,si;n2, 02,52) 



2TTi 



S2 



where the u-contour contains [—1, 1]. In this section, we prove that ^("i''^i)'("2,a2)(^g2, si) — 
r(ni,ai,si;n2, 02,52). 

Proposition 4.7. Assume {ni,ai) < {n2,a2). Then 

(r*0)"^-"i *r(s2,si) = r(ni,-l/2,si;n2,l/2,S2), 

{T*cb)"--"'{s2,si) = r(ni,l/2,si;n2,l/2,S2), 

((/)*r)"^-"i-i*0(s2,si) = r(ni,l/2,si;n2,-l/2,S2), 

(0*r)"^-"Hs2,Si) = rK,-l/2,Si;n2,-l/2,S2). 

Proof. Proceed by induction on 2n2 + 02 — (2ni + ai). First assume 2n2 + 02 — 
(2ni + fli) = 1 with ni = n2. Then 

r(ni, -1/2, si;ni, 1/2,52) 

(l-a;)-i/2(i + 2;)-i/2rfyrfa. 

X ^ ^ ^ 

X — u 

= f J(-1A-1/2)(^)j(1/2,-1/2)(^)(^_^)-1/2(^^^)-1/2^^ 



2m 



1 



j(-l/2.-l/2)(^) E W^^^'/''"'/'HOJ^'^''"'^'H^) 



r=0 



X (l-2;)-i/2(i + ^-)-i/2^^ 

= T(S2, Si). 

The second equality follows by evaluating the residues ai u = x, the third 
equality follows from Lemma 2.4 ^a), and the fourth equality follows from the 
orthogonality relations. 
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Now assume 2n2 + 02 — (2rii + ai) = 1 with ni ^ 712- Then 
r(ni, 1/2, si;ni + 1,-1/2,52) 



27rz 



(l-x)i/2(i + a;)-i/2rfyrfa; 



r j(V2-v2)(,)(j(-i/2,-i/2)(^)_,)(^_,)- 



27ri 



X (l-x)i/2(i + a;)-i/2^^ 

^(l/2,-l/2)(^^) 1 



S2-1 



2,-1/2) 



(^) 



r=0 



TT 27ri 
X (l-x)i/2(i + a;)-i/2rf2; 

= 0(s2,Sl)- 

The second equahty follows by evaluating the residues at u = x and u — I, the 
third equality follows from Lemma 2.4 'b), and the fourth equality follows from 
the orthogonality relations. 

The inductive step is proved using a similar argument with the help of Lem- 

□ 



mas 



2.4 » and (b). 



4.6 Computing the Kernel 



Let us now to compute the correlation kernel, which is given by (29 1. 

The first case is when (^2,02) < (^1,01). Then ,/,("i''ii)^("2,5^^ q. Fur- 
"2 

thermore, ni > n2, and ^'"j!!?^(si)<&"^-fc('*2) was calculated in Proposition 

fe=i 



The second case is when (ni,ai) < (7x2,02) and ni — n2- This happens only 
when ai = -1/2 and 02 = 1/2. Then -</,("i.«i)'("2,a2) ^ ^q-, which cancels 
with the single integral in ( 33 ) . 

The final case is when (ni,ai) <i (^2,02) and ni < n2- Adding Propositions 
4.6 and 4.7 we have the desired expression, plus an "extra" term: 

^(n2-ni-l)(^2^^ 



W^(°i'-V2)(g^ ) j(ai,-l/2)j(a2,-l/2) 
TT 



^(ai,-l/2)(g^) 1 

TT 27ri 



1 E{x) 
X (x - l)"i-"^(l - xY^{l + x)-^/'^dx 
1 i-E{u)~E{l) 



{n2 - ni - 1) 



(n2-"l-l)! 



(u-l) 



"2— "1—1 



X j(^--V2)(^)j(^2,-l/2)(y)(„_l)' 



(l - x)"^ (l + x)^^^^dudx 
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where the M-contour contains [—1,1] and does not contain any zeroes of E{u). 
We prove this equals zero. 

Let us evaluate the double integral. We use the identity 



E{u) - E{1) 



(n2 -ni~ 1)! 



n2 — ni — 1 



1 



E{vj)dw 



2m J (w - u)(tt; - l)"2-ni 
where the w-contour contains u and 1. The first double integral is now 



, (37) 



2 „1 



xji-'-V2)(:,)j(^..-l/2)(,) 



E{w) 
E^ 
dw 



{l~x)''^{l + x)-^^^dudx 



{W - u){w - l)"2-ni 

The u-contour has poles only at u = a;, so we get 



(l-x)''i(l + .T)-i/2du;da; 
^ (x - ■u;)(w - l)"2-»i ' ^ ' 



By ( 37 ) , the single integral equals 



(l-a;)''i(l + a;)-i/2^Wa; 
(w - x)(u) - l)"2-"i ' 



which cancels ( 38 1 . 



The proof of Theorem 4.1 is now complete. 



5 Asymptotics of the Kernel 

In this section, we analyze the large-time asymptotics of our system as 7 ^ 00. 
Figure |5] shows the result of a computer simulation of the Markov chain. Notice 
three distinct regions: one region where the particles are densely packed, another 
region where there are no particles, and an intermediate region. In section [5.1[ 
we find explicit formulas for the curves qi and (72 that separate these three 
regions. Compare Figures |5] and [6] 
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The appropriate global scaling is to take the time parameter 7 to vary pro- 
portionally to tN , while Ui and si vary proportionally to IN and dN , respec- 
tively. Assume the pairwise differences rii — nj and Si — Sj remain finite and 
constant. These limits are known as the bulk limits. In the limit N —i- 00, the 
behavior in the intermediate region is described by the incomplete beta kernel, 



which is an extension of the ubiquitous sine kernel. See Theorem 5.2 for the 
precise statement. 

There are also two other scaling limits that we consider. The first occurs 
when 7 oc tN, Ui cx IN with — Uj finite constants, and Si are finite constants. 
In other words, we are considering the large-time behavior of our point process 
at a finite distance from the wall on the left. This behavior is described by 



the discrete Jacobi kernel, which we introduce. See Theorem 5.6 for the precise 
statement. The second edge limit occurs when 7 c>c N/2, Ui cc N + rji^fN for 
some rji, and Si oc OiN^I^ for some ct;. In other words, we are zooming in at the 
point where q\ meets the y-axis in Figure |6] The behavior here is described by 
the symmetric Pearcey kernel. It is an analog of the Pearcey kernel, which has 



previously appeared in [2 [TTJ [T3 [SI HOI US]. See Theorem 5.7 for the precise 
statement. 

5.1 Bulk Limits 

Define the polynomials 

Rt^d,i{z) ^t+{2l + 2d- t)z + {21 -2d- t)z^ + tz^ 
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and 

QtAz) = l{l - '^tf - {2P + lOlt - t^)z^ + 

Let 
and 

Note that qi{t, I) and q2{t, I) only depend on t/l. They are graphed in Figure[6] 

Proposition 5.1. Assume l,d,t> 0. 

(1) Rt,d,i has two complex conjugate roots iff I ■ qi{t,l) < d < I ■ q2{t, I). 

(2) Let zq denote the nonreal root of Rt^d,i the upper-half plane, if it exists. 
Then |zo| > 1. 

(3) Let z„iax denote the largest real root of Rt^d,i- If d > I ■ q2{t,l), then 

^max ^ 1- 

(4-) Let z„iin denote the smallest real root of Rt^d,i- If d < I - qi{t,l), then 

Proof. (1) Rt^d,i{z) has nonreal roots iff its discriminant is negative. The dis- 
criminant oi Rt,d,i{z) is IQQt^d). Since Qt,i{z) = Qt,i{—z), it crosses (0,oo) at 
most two times. If t/l > 1/2, then Qt,i{Q) < and Qt,i{+'X)) = +oo, so there- 
fore Qt^i crosses (0,oo) an odd number of times. Thus Qt^i has one positive 
real root. By using the explicit formula for the roots of a quadratic polyno- 
mial, we see that this root is q2{t,l) ■ I. So in this case, Qt.i{d) is negative iff 
qi{tj) -1 = <d< q2it,l) ■ I. 

If < t/l < 1/2, then Qt,i{0) > and Qt,i{l) = lt{-16P + 13lt - 8t^) < 
and Qt,;(-l-cxD) — +00, so Qt has two positive real roots. By using the same 
formula, we see that the roots are qi{t,l) and q2{t,l). Once again, Qt^d) is 
negative iff qi{t,l) ■ I < d < q2{t, I) ■ I. 

(2) The product of the roots of Rt^d,i equals —1. So it suffices to show that 
Rt,d,i has a root in the interval (—1, 1). In fact, Rt^d,i has a root in the interval 
(-1, 0), because Rt^A-^) ^ < and Rt,d,i{0) = t > 0. 

(3) By (1), Rt,d,i has three real roots. The product of these roots is —1, and 
their sum is 1 -I- 2(d - /)/i > 1 -I- 2l{q2{t,l) - l)/t > (this follows from the 
explicit expression for q2{t,l)). We just showed in (2) that one of these roots 
is in the interval (—1,0). Therefore the sum of the other two roots is positive, 
and their product is greater than 1. This holds only if Zmax > 1- 

(4) Because qi{t,l) is positive, t/l must be less than 1/2. By (1), Rt,d,i has 
three real roots. The product of these roots is —1, and their sum is l+2{d—l) /t < 
1 + 2l{qi{t, I) — l)/t < —1. We just showed in (2) that one of these roots is in 
the interval (—1,0). Therefore the sum of the other two roots is negative, and 
their product is greater than 1. This holds only if Zmin < — 1. □ 
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Using the notation of Proposition 5.1 define zq = zo{t, d, I) to be 

{Zjnin, if d<l ■ qi{t,l) 
zo, if / • qi{t,l) < d<l- q2{t,l) 
Zmax. if d>l- q2{t,l). 

Define the function St^d,i{z) as 

St,dA^)^t^-^^+l\og(^-±^-l^-d\ogz. (39) 
Note that dSt^d,i/dz — Rt^d,i{z)/{2z^{z — 1)), so Zq is a critical point of St^d,i{z). 



Figure 6: The x-axis is d. The y-axis is I. The left curve is d = ^ • qi{\/2, 1) and 
the right curve is d— I ■ (72(1/2, 1). 




i -'-"r" I I I I I I 

1 2 3 4 ; i " 



The incomplete beta kernel B{k, I; C) is defined by (cf. [39] ) 

B{k,l;C) ^ - z)'^z-'-'dz, 

where the contour of integration crosses (0, 1) if fc > and (—00, 0) if /c < 0. 

Theorem 5.2. Let 7 and {si : I < i < k} and {rii : I < i < k} depend on N in 
such a way that 7/-/V — )■ t > and Si/N —fd>0 and rii/N —t I. Furthermore, 



AA 



assume the differences Si — Sj and rii — rij are all finite and independent on N . 
For each j, set fij = 2nj + aj — 1/2 and Sj = s-j 



a J — rij . Then 



lim pI ^{ni,ai,si; . . .■,nk,ak,Sk) 

det [B{hi - rij, Si + hi - sj - hj; zo)]'^j^^ , if qi{t, I) < d/l < q2{t, I), 

0, ' ^fd/l>q2it,l), 

[1, lfd/l<q,{t,l). 



Proof. In the expression for the correlation kernel (Theorem 4.1), let us replace 



the a;-integration over [—1, 1] by the integration over the unit circle by setting 
X — {z + z^^)/2. Let us also make the change of variable u — {v + v^^)/2. Set 
the w-contour to go outside the unit circle. The weight on [—1, 1] becomes 



dz 



(i-.r(i+.)-v^d..^rjf;i/2_,-i/2).,. 



ai 



-1/2, 



Aiz 



ai = 1/2. 



Denote the right hand side by mai{dz). Then Kjj equals 



Kjf{ni,ai,Si;n2,a2,S2) 



2ttH 



' Z-\'Z 



lyn 



j(a2 -1/2) 



V + V 



' v-\-v 



l)n2 z + Z — V — V 



—nia^ {dz)dv 



(40) 



+1 



(ni,ai)>(ri2,a2) 



I 2 



ma^{dz)\. (41) 



Lemma 5.3. 



lim (411 = ln^>h.2 



\ 27ri J\z\=const 



y I (jj ^S2+n2+a2-(si+ni+ai)-1^2^ _ 0) ^"i +"1 -2n2 -02 



(_X)ai-a2 
2ni+ai — n2 — 12 



Proof. By using (|6|, when ai = 02 = —1/2, (41 1 equals (up to 1 



2111 



\z\ = l 



Si I Si \ / S2 I — S2 



1 



("1 ,11 )!>(n2,a2)) 

dz 

z 
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Expanding the parantheses yields four terms. For the term corresponding to 
^si+s2^ deform the contour so that \z\ < 1. This will make the integral expo- 
nentially small as ^ oo. For the term corresponding to z~^^~''^, deform the 
contour so that \z\ > 1; again, this integral vanishes as — > oo. The remaining 
term is 



1 

47rz 



\z\ = l 



- 1 



1 

47ri 



\z\ = l 



- 1 



Making the substitution z ^ z ^, the second integral in (42 1 becomes 

rii — n2 



1 

2TTi 



\z\ = l 



1 



When ai = 02 = 1/2, (41 1 equals 



1 



+ _ -si-1/2 S2+1/2 _ -S2-1/2 



A/2 _ 



yl/2 _ ,-1/2 



dz. 



- 1 



(42) 
(43) 



(2-z-z-i)^. 

iiz 



J\z\ = l Z 

Making similar deformations and substitutions, we see that this expression also 
converges to (43 1. By a similar argument, when ai — —1/2 and 02 = 1/2, (41 1 
converges to 



When fli = 1/2 and a2 = —1/2, (41 1 converges to 



"1— "2 — 1 



1 



S2-si-3/2j^^l/2 _ ^-1/2-) 



1^1=1 



- 1 



dz. 



dz. 



(44) 



Using 



z+z 



1 = - — ^f; shows that the lemma holds in all cases. 



2 ^ 2 
Lemma 5.4. Assume qi{t,l) < d/l < q2{t,l). Then 



□ 



lim (|40l) -> — 



^S2-ri2-(si-ni)-l^-|^ _ 2;)^"i+'^i^^"2-a2|^^ 



( — 1)"^""^ 
2ni+ai— n2— 02 ' 



where the integration contour crosses (— c»,0). 

Proof. We only do the calculation explicitly when ai = 02 = —1/2, because the 
other cases are similar. By ([6|, 



(-1/2,-1/2) { Z + Z ^ j(-l/2,-l/2) ^ + ^ 



.,S2 I ,,-S2 



Expanding the parantheses on the right hand side yields four terms. This means 
that (40) can be written as the sum of four terms. We now proceed to evaluate 
each of these terms separately. 
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The term corresponding to z ^'i^''^ equals 



^(ai.-l/2)(^^) 

4 • 



j?^ f + 'U ^ ^ W-S2 ^ f + 1) 



1 1 - -y 



2 + Z~ 



-2 

— . (45) 



2 2 

The part of the integrand that depends on N equals exp(A^(S'(z) — S(v)\ where 
S{z) — St,d,i{z), of. (39 1. Recall that zq is a critical point of S{z). 




With the deformations as shown in Figure [7] the double integral asymptot- 



ically evaluates to zero. Since |zo| > 1 (Proposition 5.11, the contours can be 



deformed without picking up residues aX v — z ^. The residues a.t v — z are 

M^±^-l dz, (46) 



zo 



where the contour of integration crosses (— cxd,0). 

To calculate the term corresponding to z'^'^v^^, make the substitution z ^ 
z~^. Then the integrand and contour remain the same, so this term also equals 



(46) 



It remains to calculate the terms corresponding to z^^v and z "^w . 

it suffices to calculate the term correspond- 



ing to z~^i'i;~''2. Substituting v ^ v~^, the double integral again becomes (45 1, 
except now with the ?;-contour inside the unit circle. In this case, the double 
integral is asymptotically zero, because we can deform the contours as shown in 
Figure [8] without picking up any residues. 
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Figure 8: On the left is ^{S{z) — S{zo)). Shaded regions indicate > and 
white regions indicate 3? < 0. The double zero occurs at zq. 





Collecting all the terms shows that we get 



1 

27ri 



„S2-Si-1 



z + z 



ni— "2 

1 1 dz. 



□ 



Lemmas 5.3 and 5.4 prove the theorem when qi{t, I) < d/l < q2{t, I). 

Now assume d/l > q2{t,l). Once again, there are four terms in (40 1, corre- 
sponding to z^^iw*''^. First let us calculate the term corresponding to z"^^?;''^. 
This term equals ([is]). Let Zmax denote the largest real root of Rt,d,i and de- 



form the countours as shown in Figure |9[ T hen (45 1 asymptotically evaluates 
to zero. Since z^ax > 1 by Proposition |5.1| these deformations can be made 
without picking up residues at v — z~^. The residues at v — z equal 



1 



ni— 712 

1 I dz. 



where the contour goes clockwise around the origin. 

Similarly, as before, the term corresponding to z^^d^^ also equ als (|47[ ) 
the terms corresponding to z**^ 



equal zero. Thus (40) and Hit asymp- 




(47) 



and 



totically cancel out, so K'^(ni,ai,Si;nj,aj,Sj) converges to when (ni,ai) < 



(712,02). Therefore the determinant equals 0. 

When d/l < qi{t,l), the argument is similar. Let Zmin be the smallest real 
root of Rt^d,i- Make the deformations in (40 1 as shown in Figure 10 Since 
^min < ~ 1 by Proposition |5.1| these deformations can be made without picking 
up residues aX v — z~^. The integral does not pick up residues at v — z, so 



(40) converges to zero. Thus det[i4r''']^ converges to a triangular matrix. The 
diagonal entries are given by Lemma |5.3| which all evaluate to 1. Therefore the 
determinant converges 1. □ 
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5.2 Limit Shape 

Let H : M>o x 2) be the height function defined by 

i7(7, s, n) ^ \{{y, m) e C<x) : m = n,y > s}\, 

where C<g is the random point configuration of . In other words, i? (7, s, n) 
is the number of particles to the right of (s, n) at time 7. Define h to be 

h{t,d,l)^ lim ^EH{tN,[xN],[lN]) (48) 

Recall that we defined zq — ZQ{t,d,l) and S{z) — St.d.i{z) in the previous 
section. 



Proposition 5.5. The limit (48 1 exists and 



h{t,dj)^^[^]. (49) 



Proof. Note that 



EiJ(7,s,n)= ^ pj'2,(y,n), 



J/>s 
y= n+1 mod 2 

where pi is the first correlation function. Using Theorem 
nated convergence theorem, we get 



5.2 



and the domi- 



1 



lim p*^([diV],[/7V]) = -argzo 



and 



1 1 /■°° 

h{t,d,l)^ lim — Ei?(tiV, [diV], [;iV]) = — / arg zo(t, a;, Odx. 

Now take the partial derivative of the right hand side of (49 1 with respect to d. 
The result is 

/ 1 dS{zQ)\ argzo , ^^^o 

Since S'{zq) = 0, we can conclude that 

3 ( — - — I = / ^rg zo(^: a^, 0*^^ + const. 
V 27r / 27r 

Evaluating both sides at d = +oo proves that const = 0. 

□ 
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5.3 Discrete Jacobi Kernel 

For —1 < u < 1, define the discrete Jacobi kernel L(ni, ai, si, n2, 02, S2; u) on 
X X X as follows. If (nijfli) > {n2,a2), then 

L(ni,ai,Si,ri2,a2,S2;u) 



If fli) < (rt27 02), then 
L(ni,ai,si,ri2,a2,S2;") 



-1 



For ni = 712 and (ai, 02) = (—1/2, —1/2), the integral can be evaluated. Set 
V = cos~^(u). Then 

1 1 \ n ^ /sinw(si - S2) sinw(si + S2)\ 

- V 7r(si-S2) 7I"(S1 + S2) / 

When si — S2, the above expression is evaluated by L'Hopital's rule. This can 
be viewed as a discrete analog of the Bessel kernel, which arises at the hard 
edge in random matrix models, see (1.2)-(1.3) from j45] or (2.6) from ,21j . 
The discrete Jacobi kernel arises in the following limit. 

Theorem 5.6. Let 7 depend on N in such a way that j/N t > 0. Assume 
t/l > 1/2. Let si,...,Sk be fixed finite constants. Let ni,...,nk depend on 
N in such a way that Ui/N I and their differences rii — rij are fixed finite 
constants. Then 



N 



_ \ det[L{n^,ai,s.i,nj,aj,Sj;l ~ l/t)]fj^j^, if t/l > 1/2, 

' z/t/;<l/2. 

Proof. Let A{z) = tz + I ■ log(l — z). The kernel K'^(ni, ai, Si; ni, 02, S2) equals 

(50) 

+ l(n,.a,)>(n.,a.) (^^^^-l^^i^M. £^ J (^-^ V2) (^) J (a. 1/2) _ l)n.-n.(l _ ^)a,(i _^ :c)-'/^dxdu^ 

(51) 

Now deform the u contour as shown in Figure 11 With this deformation, 

^N{A{x)-A(l-l/t)) pN-'R{A{x)-A{l-l/t)) 

0, 



pAf(A(M)-A(i-;/t)) 



g7V3?(A(«)-A(l-;/*)) 
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Figure 11: On the left is ^{A{z) - A{1 - l/t)). Shaded regions indicate > 
and white regions indicate 3? < 0. The double zero occurs at 1 — l/t. 





so the integrand converges to zero. However, for t/l > 1/2, the deformations 
cause the double integral to pick up residues ai u — x. Thus, expression (50 1 
converges to 

Tl/(ai.-l/2)('„ \ /-l-i/t 

^ / Ji^i'-i/2)(x)ji^^-i/2)(x)(x-l)"^-"^(l-a;)'^Hl+x)-i/2rfx. 

(52) 

Adding (51) to (52 1 shows that K'^ converges to the discrete Jacobi kernel. 

U t/l < 1/2, then the double integral does not pick up residues a,t u = x. 
Thus det[_ftr'''] converges to a triangular matrix. The diagonal entries are given 
by (51), which all evaluate to 1. Therefore the determinant equals 1. □ 



5.4 Symmetric Pearcey Kernel 



Define a kernel /C on x M as follows. In the expressions below, the u-contour 
is integrated on rays from ooe*'^/^ to to ooe~"/^. If t?i < 772, then 

/C(cri,77i,(T2,'72) = ^ / / exp(-77ia;^ + 772U^ + - a;"*) 



u 

X cos{aix) cos{a2u) — -dxdu. (53) 

— 



If 772 < I then 



u-i ^ I {(Tl + <J2f , (cti-CT2)2 
/C(cri,77l,(T2,7?2) = , , °^PT7 V +CXP 



2^7r(77i - 772) V 4(772-771) 4(7^2-771). 
exp(— 7710;^ + ri2U^ + m"* — x"^) cos(cria;) co&{(J2u)^ Tjdxdu. 

— X'^ 

(54) 
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This kernel arises as follows. Let pj,'^ denote the correlation function of 
Vy^ A' which is the pushforward of 7^ J under A. See Appendix [Aj 

Theorem 5.7. For \ < i < k, let Si depend on N in such a way that s,/N^/^ 
2"^/*cri > as iV — > oo. Let 7 depend on N in such a way that j/N — s- 1/2 
as N ^ 00. Let Ui depend on N in such a way that [rii — N)/^/N —> 2^^^^rii. 
Then 

lirn^ i -^jrj^ j p1'^{ni,ai,si; . . . ; n^, afc, s^) = det [/C (ctj, 77,^, Cj, ?7J)]l<^J<fc■ 
Proo/. From Corollary |4.2| the left hand side of the above equation equals 
det[(iVi/4/25/4)^^(„^^a„s„nj,a„5j)]^_^-^iand {N^'y2''/^)Kl{ni,ai,si,n2,a2, 
equals 



^'-'^^'^'^'^^'''d.du (55) 
(m - u-x ^ ^ 



- 1)"!-"^ (1 - (1 + x)-'''dx J . (56) 



Deform the contours as shown in Figure 11 with the double critical point 
at —1. Then, asymptotically, nonvanishing contributions to (55 1 and (56 1 come 
from near —1. This justifies the substitutions x' — N^^^{x + 1) and u' = 
Afi/2(u+l). For large N, x' is integrated from to cx3 and u' is integrated from 
ioo to ~ioo. There are also the following asymptotic relations: 

{l-x)'''il + x)-^^^ w 2'''N^/\x'y^/^, (57) 
dxdu _ du'dx' 
u — X u' — x' 

Let us show that if sjN'^l^ then 

(_l).j(±i/2,-i/2)(^) -> cos(2-3/4^(a;')i/2) (59) 

For fli = -1/2, 



Hence, 



j(-l/2,-l/2)(^) = C0S(S6'), X = COS (9. 

(-l)^Jiri/2'-i/2)(a;) ^ (-1)^ • cos(s • cos-i(-l + N-^/'^x')) 
(-1)'' • cos(s7r - 2-3/V(x')i/^ + o(l)) 

(-l)^(cos(s^) cos(2-3/4a(x')'/') + sin(s^) sin(2-3/V(:E')'^')) 
cos(2-3/4a(xy/2) 
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Similarly, for a~ 1/2, 



sin((s + 1/2)6) 



sin(6l/2) 

sin((s + l/2)7r - 2-3/V(2;')^/^ + o(l)) 



V(2 + o(l))/2 

« (-l)"(sin((s + i)^) cos(2-3/V(x')'/') - cos((s + \)tt) s\n{2-^'^a{x')^''')) 
= cos(2-3/V(a;')'/') 

For fli = ±1/2 and s > 0, 

= — (60) 



Let A{z) = z/2 + log(l - z). We have 



- A{~1) = --(z + 1)2 + 0((z+ 1)3), 



and 



(_2)"^~"i exp(7(x-M)) 



(a; - l)"i gAr(A(a:)--A(-l))+2-i''^,,iV]V(-log2+log(l-a;)) 
(U-1)"2 gAf(A(«)-A(-l))+2-i/2,,2y]v(-log2+log(l-a:)) 

g-(2;')^/8-»?l2;72'/^ 



-(«')V8-»?2«72^/^ ■ 



(61) 



The kernel can be multiplied by the conjugating factor (— 1)**^ '''^[—2)'^'^ ""^ 
without changing the determinant. Combining (551 — (61) shows that 



(_X)si-s2(_2)'- 



^(551 



_/Vl/** (_l)si+S2('_2)"2-ni ^(ai,-l/2) 



25/4 



2m 



|u| = l J-1 



j(--V2)(^)j(.2,-l/2)(^) 



■u Si 



(x-l)^(l-a;)'^i(l + a;)-i/2 
X ^ TTTT ^ ^ — ^ dxdu 



{u-l) 



N 



-zoo /'OO 



exp 



21/4 



™'-l/2 

X COS {2'^'^a2u'^'^)-, -dx'du' 

w — x' 



(-^^1^' + ^2^^^' + ^(-" - ^")) cos(2-3/Vx'i/2) 



= ([531. 

In the last equality, the substitutions u' — 2?l'^u^ and x' = 2'^/2i2 were needed. 
We need one final additional calculation: 

(_2)"-»i(a:-l)"i-"^ «(l-i7V-l/V)(''l-''^)2-V2^/]v^^^p(_^(^^„^^)^,)^ 
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which shows that 



(-l)"i-"'^(-2)"2-"i(|56l 

2I/4 roo 



exp 



'-ni>V2 



(cti+CT2)^ (cri-cr2)^ 

/ . exp — + exp — 

2V7r(77i - 772) V 4(?72-?7i) 4(772-771) 



Therefore (-1)^1-"^ (-2)"^-"! (d55| + (561) f3 



□ 



A Generalities on Random Point Processes. 

Let X be a locally compact separable topological space. A point configuration 
X in X is a locally finite collection of points of the space X. For our purposes 
it suffices to assume that the points of X are always pairwise distinct. Denote 
by Conf (X) the set of all point configurations in X. 

A relatively compact Borel subset ^ C X is called a window. For a window 
A and X G Conf(X), set Na{X) = [A n X| (number of points of X in the 
window). Thus, Na is a function on Conf(X). Conf(X) is equipped with the 
Borel structure generated by functions Na for all windows A. 

A random point process on X is a probability measure on Conf(X). One 
often uses the term particles for the elements of a random point configuration. 

Given a random point process on X, one can usually define a sequence 
where pn is a symmetric measure on X" called the rtth correlation 
measure. Under mild conditions on the point process, the correlation measures 
exist and determine the process uniquely. 

The correlation measures are characterized by the following property: For 
any tt, > 1 and a compactly supported bounded Borel function / on X" one has 



fPn = 




where ( • ) denotes averaging with respect to our point process, and the sum 
on the right is taken over all n-tuples of pairwise distinct points of the random 
point configuration X. 

Often one has a natural measure /i on X (called reference measure) such that 
the correlation measures have densities with respect to /^^", n = 1,2, ... . Then 
the density of p„ is called the nth correlation function and it is usually denoted 
by the same symbol p„. 

The first correlation function pi is often called the density function as it 
measures the average density of particles. 

For point processes on a finite or countable discrete space X it is natural to 
choose the counting measure as the reference measure /x, and then there is a 
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simpler way to define the correlation functions: For any n = 1,2, .. . and any 
pairwise distinct xi, . . . , a;„ G X, 

Pn{xi, ... ,x„) = Prob{X e Conf (X) \ X D {xi, . . . , x„}}. 

If X is discrete, a random point process on X is always uniquely determined 
by its correlation functions. 

The reader can find more information on random point processes in |17| . 

A point process on X is called determinantal if there exists a function K{x, y) 
on X X X such that the correlation functions (with respect to some reference 
measure) are given by the determinantal formula 

Pn{xi, . . . ,x„) = det[K{x.„Xj)X^,^^^ 

for all n = 1, 2, . . . . The function K is called the correlation kernel. 

Note that the correlation kernel is not defined uniquely: K{x, y) and j^K{x, y) 
define the same correlation functions for an arbitrary nonzero function / on X. 

Assume that X is discrete. Define a map A by 

A : Conf(X) ^ Conf(X), X ^ X\X. 

Given a point process V on X, its pushforward under A is also a point process on 
X; denote it by Va- The map A is often referred to as particle-hole involution., 
because the particles of Va are located exactly at those points of X where there 
are no particles of V . With this notation, we have the following proposition. 

Proposition A. // V is a determinantal point process with correlation kernel 
K(x,y), then Va is also a determinantal point process with correlation kernel 

KA{x,y) = S.^^y - K{x,y). 

The proof is an application of the inclusion-exclusion principle, see Proposi- 
tion A.8 of [13]. 
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